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We propose an improved scheme of perturbation theory based on our exact solution [An Min 
Wang, quant-ph/0602055] in general quantum systems independent of time. Our elementary start- 
point is to introduce the perturbing parameter as late as possible. Our main skills are Hamiltonian 
redivision so as to overcome a flaw of the usual perturbation theory, and the perturbing Hamiltonian 
matrix product decomposition in order to separate the contraction and anti-contraction terms. Our 
. . . calculational technology is the limit process for eliminating apparent divergences. Our central idea 

' is "dynamical rearrangement and summation" for the sake of the partial contributions from the high 

' order even all order approximations absorbed in our perturbed solution. Consequently, we obtain 

, the improved forms of the zeroth, first, second and third order perturbed solutions absorbing the 

partial contributions from the high order even all order approximations of perturbation. Then we 
^ ' deduce the improved transition probability. In special, we propose the revised Fermi's golden rule. 

Q , Moreover, we apply our scheme to obtain the improved forms of perturbed energy and perturbed 

• state. In addition, we study an easy understanding example of two-state system to illustrate our 

' scheme and show its advantages. All of this implies the physical reasons and evidences why our 

CNj , improved scheme of perturbation theory are actually calculable, operationally efhcient, conclusively 

04 ■ more accurate. Our improved scheme is the further development and interesting application of our 

exact solution, and it has been successfully used to study on open system dynamics [An Min Wang, 
CN ' quant-ph/0601051]. 
> ' 

' PACS numbers: 03.65.-w, 04.25.-g, 03.65. Ca 
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I. INTRODUCTION 



The known perturbation theory P, Q is an extremely important tool for describing real quantum systems, as it 
turns out to be very difficult to find exact solutions to the Schrdinger equation for Hamiltonians of even moderate 
^ complexity. Recently, we see the dawn to overcome this difficulty because we obtained the exact solution of the 

Schroding equation pj] in general quantum systems independent of times Q. However, this does not mean that the 
^ , perturbation theory is unnecessary because our exact solution is still an infinite power series of perturbation. Our 
solution is called "exact one" in the sense including all order approximations of perturbation. In practice, if we do 
not intend to apply our exact solution to investigations of the formal theory of quantum mechanics, we often need to 
cut off our exact solution series to some given order approximation in the calculations of concrete problems. Perhaps, 
one argues that our exact solution will back to the usual perturbation theory, and it is, at most, an explicit form 
that can bring out the efficiency amelioration. Nevertheless, the case is not so. Such a view, in fact, ignores the 
significance of the general term in an infinite series, and forgets the technologies to deal with an infinite series in the 
present mathematics and physics. From our point of view, since the general term is known, we can systematically 
and reasonably absorb the partial contributions from some high order even all order approximations to the lower 
order approximations just like one has done in quantum field theory via summation over a series of different order 
but similar feature Feynman figures. In this paper, based on such a method we develop our "dynamical arrangement 
and summation" idea, and then propose an improved scheme of perturbation theory via introducing several useful 
skills and methods. 

It is very interesting that we find a flaw in the usual perturbation theory, that is, the perturbing parameter is 
introduced too early so that the contributions from the high order even all order approximations of the diagonal and 
off-diagonal elements of the perturbing Hamiltonian matrix are, respectively, inappropriately dropped and prematurely 
cut off. For some systems, the influences on the calculation precision from this flaw can be not neglectable with the 
evolution time increasing. This motivates us to set our start-point to introduce the perturbing parameter as late as 
possible in order to guarantee the generality and precision. It is natural from a mathematics view if we think the 
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perturbing parameter in a general perturbation theory is a formal multiplier. Based on this start point, we propose 
Hamiltonian redivision skill and further methods so as to overcome the above flaw in the usual perturbation theory, 
viz. the Hamiltonian redivision makes the contributions from all order approximations of the diagonal elements of 
the perturbing Hamiltonian matrix can be absorbed in our improved form of perturbed solution. Hence, this skill 
advances the calculation precision in theory, extends the application range for the perturbation theory and can remove 
degeneracies in some systems. 

Since our exact solution series has apparent divergences, we provide the methods of perturbing Hamiltonian matrix 
product decomposition in order to separate the contraction terms with apparent divergences and anti-contraction 
terms without apparent divergences. Here, "apparent" refers to an unture thing, that is, the apparent divergences are 
not real singularities and they can be eliminated by mathematical and/or physical methods, while the "perturbing 
Hamiltonian matrix" refers to the representation matrix of the perturbing Hamiltonian in the unperturbed Hamilto- 
nian representation. Then, by the limit process we can eliminate these apparent divergences in the contraction terms. 
Furthermore we apply "dynamical rearrangement and summation" idea for the sake of the partial contributions from 
the high order even all order approximations absorbed in our perturbed solution. In terms of these useful ideas, skills 
and methods we build an improved scheme of perturbation theory. Without any doubt, they are given definitely 
dependent on our exact solution. In fact, our exact solution inherits the distinguished feature in a c- number function 
form just like the Feynman 0] path integral expression and keeps the advantage in Dyson series jsj that is a power 
series of perturbation. At the same time, our exact solution is so explicit that when applying it to a concrete quantum 
system, all we need to do is only the calculations of perturbing Hamiltonian matrix and the limitations of primary 
functions. 

As well known, a key idea of the existed perturbation theory to research the time evolution of system is to split the 
Hamiltonian of system into two parts, that is 

H = Ho + Hi, (1) 

where the eigenvalue problem of so called unperturbed Hamiltonian Hq is solvable, and so-call perturbing Hamiltonian 
Hi is the rest part of the Hamiltonian. In other words, this splitting is chosen in such a manner that the solutions of 
Hq are known as 

Fol*^) --B^l*^): (2) 

where 1$''') is the eigenvector of Hq and is the corresponding eigenvalue. Whole \^^), in which 7 takes over all 
possible values, form a representation of the unperturbed Hamiltonian. It must be emphasized that the principle of 
Hamiltonian split is not just the best solvability mentioned above in more general cases. If there are degeneracies, 
the Hamiltonian split is also restricted by the condition that the degeneracies can be completely removed via the 
usual diagonalization procedure of the degenerate subspaces and the Hamiltonian redivision proposed in this paper, 
or specially, if the remained degeneracies are allowed, it requires that the off-diagonal element of the perturbing 
Hamiltonian matrix between any two degenerate levels are always vanishing so as to let our improved scheme of 
perturbation theory work well. As an example, it has been discussed in our serial study In addition, if the 
cut-off approximation of perturbation is necessary, it requires that for our improved scheme of perturbation theory, 
the off-diagonal elements of Hi matrix is small enough compared with the diagonal element oi H ~ Hq + Hi matrix 
in the unperturbed representation. 

Nevertheless, there are some known shortcomings in the existed perturbation theory, for example, when Hi is 
not so small compared with Hq that the high order approximations should be considered, and/or when the partial 
contributions from the higher order approximations become relatively important to the studied problems, and/or the 
evolution time is long enough, the usual perturbation theory might be difficult to calculate to an enough precision 
in an effective way, even not feasible practically since the lower approximation might break the physical symmetries 
and/or constraints. In order to overcome these difficulties and problems, we recently study and obtain the exact 
solution in general quantum systems via explicitly expressing the time evolution operator as a c-number function 
and a power series of perturbation including all order approximations [3||. In this paper, our purpose is to build 
an improved scheme of perturbation theory based on our exact solution '3\ so that the physical problems are more 
accurately and effectively calculated. For simplicity, we focus on the cases of Schrodinger dynamics jlj] . It is direct to 
extend our improved scheme to the cases of the von Neumann dynamics Q • 

Just well-known, quantum dynamics and its perturbation theory have been sufficiently studied and have many 
successful applications. Many famous physicists created their nice formulism and obtained some marvelous results. 
An attempt to improve its part content or increase some new content as well as some new methods must be very 
difScult in their realizations. However, our endeavors have obtained their returns, for examples, our exact solution 
0, perturbation theory and open system dynamics @ in general quantum systems independent of time. 

In this paper, we start from proposing our ideas, skills and methods. We expressly obtain the improved forms of the 
zeroth, first, second and third order approximations of perturbed solution absorbing partial contributions from the 



high order even all order approximations, finding the improved transition probability, specially, the revised Fermi's 
golden rule, and providing an operational scheme to calculate the perturbed energy and perturbed state. Furthermore, 
by studying a concrete example of two state system, we illustrate clearly that our solution is more efficient and more 
accurate than the usual perturbative method. In short, our exact solution and perturbative scheme are formally 
explicit, actually calculable, operationally efficient, conclusively more accurate (to the needed precision). 

This paper is organized as the following: in Sec. [ill we find a flaw of the usual perturbation theory and introduce 
Hamiltonian redivision to overcome it. Then, we propose the skill of the perturbing Hamiltonian matrix product 
decomposition in order to separate the contraction terms with apparent divergences and anti-contraction terms without 
apparent divergences. By the limit process we can eliminate these apparent divergences. More importantly, we use 
the "dynamical rearrangement and summation" idea so that the partial contributions from the high order even all 
order approximations are absorbed in our perturbed solution and the above flaw is further overcome; in Sec. IIIII we 
obtain the improved forms of the zeroth, first, second and third order perturbed solutions of dynamics absorbing 
partial contributions from the high order even all order approximations; in Sec. II VI we deduce the improved transition 
probability, specially, the revised Fermi's golden rule. In Sec. |V] we provide a scheme to calculate the perturbed 
energy and the perturbed state; in Sec. IVII we study an example of two state system in order to concretely illustrate 
our solution to be more effective and more accurate than the usual method; in Sec. I VIII we summarize our conclusions 
and give some discussions. Finally, we write an appendix as well as a supplementary where some expressions are 
calculated in order to derive out the improved forms of perturbed solutions. 



II. SKILLS AND METHODS IN THE IMPROVED SCHEME OF PERTURBATION THEORY 



In our recent work [3|, by splitting a Hamiltonian into two parts, using the solvability of eigenvalue problem of 
one part of the Hamiltonian, proving an useful identity and deducing an expansion formula of operator binomials 
power, we obtain an explicit and general form of the time evolution operator in the representation of solvable part 
(unperturbed part) of the Hamiltonian. Then wc find out an exact solution of the Schrodinger equation in general 
quantum systems independent of time 
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and all ~ (<I>'''j form so-called "perturbing Hamiltonian matrix", that is, the representation matrix 

of the perturbing Hamiltonian in the unperturbed Hamiltonian representation. While 
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here 2 < i < I. 

It is clear that there are apparent divergences in the above solution. For example 
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when = E^i (which can appear in the summation or degeneracy cases), it is oo — oo = —iHl^ t. As pointed out 
in our paper ^3|], we need to imderstand AJ^ (t) in the sense of hmitation. Moreover, in practice, we should present 
how to calculate their limitation in order to eliminate the apparent divergences. 

Now, the key problems are how and when to introduce the cut-off approximation in order to obtain the perturbed 
solution. For studying and solving them, we first need to propose some skills and methods in this section. These skills 
and methods profit from the fact that the general term is clearly known and explicitly expressed in our exact solution. 
By using them we can derive out the improved forms of perturbed solution absorbing the partial contributions from 
the high order even all order approximations of perturbation. It will be seen that all the steps are well-regulated 
and only calculational technology is to find the limitation of primary functions. In other words, our exact solution 
and perturbation theory are easily calculative and operational, and they have better precision and higher efficiency. 
Frankly speaking, before we know our exact solution, we are puzzled by too many irregular terms and very trouble 
dependence on previous calculation steps. Moreover, we are often anxious about the precision of the results in such 
some calculations because those terms proportional to t"e~''^"'i* (a = 1,2, ••■) in the high order approximations 
might not be ignorable with time increasing. Considering the contributions of these terms can obviously improve the 
precision. However, the known perturbation theory does not give the general term, considering this task to absorb 
reasonably the high order approximations is impossible. 

Since our exact solution has given the explicit form of any order approximation, that is a general term of an 
arbitrary order perturbed solution, and their forms are simply the summations of a power series of the perturbing 
Hamiltonian. Just enlightened by this general term of arbitrary order perturbed solution, we use two skills and 
"dynamical rearrangement and summation" method to build an improved scheme of perturbation theory, which are 
respectively expressed in the following two subsections. 



A. Hamiltonian Redivision 



The first skill starts from the decomposition of the perturbing Hamiltonian matrix, that is the matrix elements of 
Hi in the representation of Hq, into diagonal part and off-diagonal part: 



737j + l 



„7j73 + l 



(11) 



SO as to distinguish them because the diagonal and off-diagonal elements can be dealed with in the different way. In 
addition, it makes the concrete expression of a given order approximation can be easily calculated. Note that hj^ has 
been chosen as its diagonal elements and then gl^''^'^^ has been set as its off-diagonal elements: 



„7j7j + l 
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(1 Sj-j-^-^). 

As examples, for the first order approximation, it is easy to calculate that 
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Note that here and after we use the symbol AT* denoting the contribution from the ith order approximation, which 
is defined by ©, while its argument indicates the product form of perturbing Hamiltonian matrix. However, for the 
second order approximation, since 
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we need to calculate the mixed product of diagonal and off-diagonal elements of perturbing Hamiltonian matrix. 
Obviously, we have 
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In the usual time-dependent perturbation theory, the zeroth order approximation of time evolution of quantum state 
keeps its original form 
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where we have set the initial state as \^^) for simplicity. By using our solution, we easily calculate out the contributions 
of all order approximations from the product of completely diagonal elements h of the perturbing Hamiltonian matrix 
to this zeroth order approximation 
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Therefore, we can absorb the contributions of all order approximation parts from the product of completely diagonal 
elements h of the perturbing Hamiltonian matrix into this zeroth order approximation to obtain 

^''W(t)^ =e-'(^-^+''?)*|i>T). (22) 

Similarly, by calculation, we can deduce that up to the second approximation, the perturbed solution has the following 
form 
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However, for the higher order approximation, the corresponding calculation is heavy. In fact, it is unnecessary to 
calculate the contributions from those terms with the diagonal elements of Hi since introducing the following skill. 
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This is a reason why we omit the relevant calculation details. Here we mention it only for verification of the correctness 
of our exact solution in this way. 

The results ([22|) and ([23|) are not surprised because of the fact that the Hamiltonian is re-divisible. Actually, we 
can furthermore use a trick of redivision of the Hamiltonian so that the new Hq contains the diagonal part of Hi, 
that is, 

= Ho+Y,hJ\^-'){<^'>\, (24) 

7 

H[ = Hi-J2hJ\^^){<S>'f\=Y,gr'm{'^n (25) 

7 7.7' 

In other words, without loss of generality, we always can choose that H[ has only the off-diagonal elements in the Hq 
(or Hq) representation and 

jT^I$^) = (£;^ + /i7)l<i>^) = £;;i<i>^). (26) 

It is clear that this redivision does not change the representation of the unperturbed Hamiltonian, but can change the 
corresponding eigenvalues. In spite that our skill is so simple, it seems not be sufficiently transpired and understood 
from the fact that the recent some textbooks of quantum mechanics still remain the contributions from the diagonal 
elements of the perturbing Hamiltonian matrix in the expression of the second order perturbed state. It is clear 
that the directly cut-off approximation in the usual perturbation theory drops the contributions from all higher order 
approximations of the diagonal element of the perturbing Hamiltonian matrix. From our point of view, the usual 
perturbation theory introduces the perturbing parameter too early so that this flaw is resulted in. 
If there is degeneracy, our notation has to be changed as 

^Ji ~^ ^Ha-ti =-^7ii (27) 
'^7i73 * '^7i7j '^i^.a^^. I (28) 
^7i73 ^ Vjilj ~^ ^JiljVa^.a^ ■ (29) 



Thus, we can find 
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This seems to bring some complications. However, we can use the trick in the usual degenerate perturbation theory, 
that is, we are free to choose our base set of unperturbed kets |$^°^) in such a way that that H\ is diagonalized in 
the corresponding degenerate subspaces. In other words, we should find the linear combinations of the degenerate 
unperturbed kets to re-span the zero-order eigen subspace of Hq so that 



($7a.|^^|$7^) ^ gl^^''^''^ ^dl^^6,_^„^. 



(32) 



(If there is still the same values among all of d-ya, this procedure can be repeated in general.) This means that 
^70^70^' _ Then, we use our redivision skill again, that is 



7^-D, 7 a-y 



(33) 
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where 13 is a set of all degenerate subspace-indexes. Thus, the last term in Eq. ([?T|) vanishes, 
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In fact, under the preconditions of Hi is diagonal in the degenerate subspaces, we can directly do replacement 

gr ^ gr'""'^^' V,n, (36) 

from the non-degenerate case to the degenerate case. For simplicity, we always assume that Hi has been diagonalized 
in the degenerate subspaces from now on. 

It must be emphasized that the Hamiltonian rcdivision skill leads to the fact that the new perturbed solution can 
be obtained by the replacement 

E^, E-y, + hi' (37) 

in the non-degenerate perturbed solution and its conclusions. With degeneracy present, if our method is to work well, 
the degeneracy should be completely removed in the diagonalization procedures of the degenerate subspaces and the 
Hamiltonian rcdivision, that is, for any given degenerate subspace, dl"" ^ dj^ if a ^ 6. In other words, E'^^ ^ E'^j^ 
ii a b. This means that all of eigenvalues of Hq are no longer the same, so we can back to the non-degenerate 
cases. Or specially, if we allow the remained degeneracies, the off-diagonal element of the perturbing Hamiltonian 
matrix between any two degenerate levels are always vanishing. This implies that there is no extra contribution from 
the degeneracies in the any more than the zeroth order approximations. It is important to remember these facts. 
However, how must we proceed if the degeneracies are not completely removed by the usual diagonalization procedure 
and our Hamiltonian rcdivision as well as the special cases with the remained degeneracies stated above are not valid. 
It is known to be a challenge in the usual perturbation theory. Although our exact solution can apply to such a 
kind of cases, but the form of perturbed solution will get complicated because more apparent divergences need to 
be eliminated and then some new terms proportional to the power of evolution time will appear in general. We will 
study this problem in the near future. Based on the above reasons, we do not consider the degenerate case from now 
on. 

From the statement above, we have seen that there are two equivalent ways to obtain the same perturbed solution 
and its conclusions. One of them is to redefine the energy level i?-y. as E!y. (or £"".), think E'^. (or i?".) to be explicitly 
independent on the perturbing parameter from a redefined view, and then use the method in the usual perturbation 
theory to obtain the result from the redivided H[ (or Hi). The other way is to directly derive out the perturbed 
solution from the original Hamiltonian by using the standard procedure, but the rearrangement and summation are 
carried out just like what we have done above. From our point of view, this is because the perturbing parameter is 
only a formal multiplier in mathematics and it can be introduced after redefining E!y. . It is natural although this 
problem seems not be noticed for a long time. The first skill, that is, the Hamiltonian rcdivision skill will be again 
applied to our scheme to obtain improved forms of perturbed energy and perturbed state in Sec. IIIII 

The Hamiltonian rcdivision not only overcomes the flaw of the usual perturbation theory, but also has three obvious 
advantages. Firstly, it advances the calculation precision of perturbation theory because it makes the contributions 
from all order approximations of the diagonal elements of the perturbing Hamiltonian matrix are naturally included. 
Secondly, it extends the applicable range of perturbation theory based on the same reason since the diagonal elements 
of the perturbing Hamiltonian no longer is needed to be smaller. Lastly, it can be used to remove the degeneracies, 
which is important for the perturbation theory. 

For simplicity, in the following, we omit the ' (or ") in Hq, Hi as well as E^, and always let Hi have only its 
off-diagonal part and let Hq have no degeneracy unless particular claiming. 



B. Perturbed Hamiltonian matrix product decomposition and apparent divergence elimination 

In this subsection, we present the second important skill enlightened by our exact solution, that is, the perturbing 
Hamiltonian matrix product decomposition, which is a technology to separate the contraction terms with apparent 
divergences and anti-contraction terms without apparent divergences, and then we can eliminate these apparent diver- 
gences by the limit process. More importantly, we can propose so-called "dynamical rearrangement and summation" 
idea in order to absorb the partial contributions from the high order even all order approximation of perturbing 
Hamiltonian into the lower order terms of our perturbation theory. It is a key method in our improved scheme of 
perturbation theory. 

Let us start with the second order approximation. Since we have taken h'I''''^^ only with the off-diagonal part 

^7jTj+i^ the contribution from the second order approximation of the perturbing Hamiltonian is only A^'^ (gg) in 

eg. p^ . However, we find that in the above expression of AJ^ (gg), the apparent divergence has not been completely 
eliminated or the limitation has not been completely found out because we have not excluded the case Ey — Ej/ (or 
7 = 7'). This problem can be fixed by introducing a perturbing Hamiltonian matrix product decomposition 

„7l72 „72 73 _ „7l 72^72 73 1 _l_ ^7l72 „72 73„ (oo\ 
gi i/l — i/l i/l "7173 "T i/l i/l '/7l73' \'~'°) 
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where 77-^^- 
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Thus, the contribution from the second order approximation is made of two terms, one 



so-called contraction term with the 6 function factor and the other so-called anti-contraction term with the rj function 
factor. Obviously, the contraction term has the apparent divergence and anti-contraction term has no the apparent 
divergence. Hence, in order to eliminate the apparent divergence in the contraction term, we only need to find its lim- 
itation. It must be emphasized that we only consider the non-degenerate case here and after for simplification. When 
the degeneration happens, two indexes with the same main energy level number will not have the anti-contraction. 

In terms of the above skill, we find that the contribution from the second order approximation is made of the 
corresponding contraction- and anti-contraction- terms 
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The above method can be extended to the higher order approximation by introducing a skill of perturbing Hamil- 
tonian matrix product decomposition, or simply called it g'-product decomposition when the perturbing Hamiltonian 
matrix is off-diagonal. For a sequential product of off-diagonal elements g with the form HfeLi 57'°''^'°^^ ("^ > 2), we 
define its (m — l)th decomposition by 



n 9?'-^' = n ar^^' '^7i7.+i + n ar-^' ^7..™+- 



(42) 



fe=i 



U=l 



When we calculate the contributions from the nth order approximation, we will first carry out (n — 1) the first 
decompositions, that is 



„7fc7fc+l 
i/l 



Kk=l 



n-1 

Y\ ('^7fc7(=+2 + %fc 7(1+2 

Lfe=i 



Obviously, from the fact that Hi is usually taken as Hermit one, it follows that 
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(44) 



When the contribution from a given order approximation is considered, the summation over one of two subscripts will 
lead to the contraction of ^-production. More generally, for the contraction of even number 51-production 



I m—l \ /m— 1 \ 

n -"^l^ n '^7fc7fc + 2 j '^7l7^7m+l7' ~ Ifl'l''' I I ^7fc7)!+2 1 ^7l7^7m+l7' ^77' ' 
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and for the contraction of odd number g'-production, 

I m — l 



gi^ ' I '^7l7'^7m+l7' 



gV 
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where ^7i7^7„+i7' is a factor appearing in the expression of our solution. 



(46) 
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Then, we consider, in turn, all possible the second decomposition, the third decomposition, and up to the (n— l)th 
decomposition. It must be emphasized that after calculating the contributions from the terms of lower decompositions, 
some of terms in the higher decompositions may be trivial because there are some symmetric and complementary 
symmetric indexes in the corresponding results, that is, the products of these results and the given or rf-y^^^, 

are zero. In other words, such some higher decompositions do not need to be considered. As an example, let us 
analyze the contribution from the; third order approximation. It is clear that the first decomposition of a sequential 
production of three off-diagonal elements becomes 



„7172„7273 _7374 _ „7l72„7273„7374 r r , „7l72 „7273 ^7374 r „ 

— ill .Vl i/l "71 73 "72 74 ^ .VI yi i/l "7l73'/7274 



9i ' 9i "91 



I 7172 7273 7374 c I 7172 7273 7374 

<yi 91 9l '/7i73"7274 I i/l 9l 9l '/7i73'/7274- 

Thus, the related contribution is just divided into 4 terms 

{999) = Af' {999; cc) + Ap {ggg; cn) + Ap' {ggg; nc) + Ap {ggg, nn). 



(47) 



(48) 



In fact, by calculating we know that the second decomposition of the former three terms do not need to be considered, 
only the second decomposition of the last term is nontrivial. This means that 



(999; nn) = Ap {ggg; nn, c) + Ap {ggg; nn, n). 



(49) 



where we have added in the definition of A^'^ {ggg;nn,c), and 77-Y173 in the definition of Ap {ggg;nn,n). Ob- 

viously, in the practical process, this feature largely simplifies the calculations. It is easy to see that the number of 
all of terms with contractions and anti-contractions is 5. For convenience and clearness, wc call the contributions 
from the different terms in the decomposition of (7-product as the contractions and anti-contractions of (7-product. Of 
course, the contraction and anti-contraction refer to the meaning after summation(s) over the subscript(s) in general. 
Moreover, here and after, we drop the argument gg • ■ ■ g in the ith order approximation Ai since its meaning has 
been indicated by i after the Hamiltonian is redivided. For example, the explicit expressions of all contraction- and 
anti-contraction terms in the third order approximation A3 can be calculated as follows: 



AP (cc) 



E 



7i.--- >74 Li=l 



d,(£;[7,3]) 
2e-'^v* 



„7.773 + l 



9i 



) <57i7'^747' 



\k=l 



3 H 3 + (-i*) 2 + 2 



9r 



„77 

9i ) 



(50) 



AP {cn) 



7l.--- >74 



E 
E 



E(-ir^ 



71 



E'yt 



„7j73 + l 



9i 



i=i 



'^71 73 ^72 74 ^71 7 ^71 + 1 " 



g— ii?^^ t 



{Ej — Eji) {E.y — Eji)"^ {E^ — Ejj^)'^ {E.y — E^i) {E.y — £/-yj)^ (-^71 ~ -^7') 



g-iiJyt 



i^E^ E'y/'j (-£'71 -^7') 



+ (-it) 



^7 71 / V 7 7 

Q — iE-yt 

{E^ — -^71) (^7 — -^7') 



\9rf9r'v,.Y 



(51) 



AP'{nc) = J2 

71, ■■ 

= E 



E(-ir^ 



di{E[j,3]) 



^771 73 ^72 74 ^71 7 ^7i + 1 1 



71 



{Ej — ^71) (-^7 — -^7')'^ (-^7 — -^71) (-^71 — -^7')^ i-^-f — ^71) (-^7 — -^7')^ 



(-^7 - -^71) (-^71 - -^7')' 



+ (-it) 



{E^ - Ey) {E^^ - Ey) 



77 

gP 



9r 



^771 



(52) 
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Ap {nn, c) 



E 

7ir-- 174 L 

E 

7172 



i=l 



di(i?[7,3]) 



n 



„7j73 + l 

yi 



^7l7 "^747' ^71 73 '?72 74 "^77' 



+ 



+ 



{E^ - Ejj'^ {E^^ - E^^ ) {E^ - -2^72 )^ (-^71 ~ -^72 ) {E-y - E^if {E^i - E^2 ) 



(-£^7 ~ -^72) (-^71 ~ -^72) 



i-it) 



Q — iEyt 



{Ej - E^^) {Ej^ - E^J 



^771^7172^727 ^^^j^^,^ (53) 



A3'*' {nn,n) 



= E 
71, •• 

= E 



E(-ir' 



,i=l 



rfi(i^[7,3]) 



'^717 "^74 7' ^71 73 ^72 74 ^77 ' 



(-£^7 ~ -^71) (-^7 ^ -^72) (-^7 ~ Ey) {E^ - Ej^) {E^^ - Ej^) [E^^ - Ey) 



+ 



{E^ - Ej^) (-E71 - £^72) (£^72 - Ey) {E^ - Ey) {Ej^ - Ey) {E^2 - Ey) 



^ 91'^^ 91^'^^ ^772 ^i 7' ^7' 



(54) 



In the above calculations, the uesed technologies mainly to find the limitation, dummy index changing and summation, 
as well as the replacement S'i'^'%i7, = gj''^^ since gj^'^' has been off-diagonal. 

It must be emphasized that, in our notation, AJ'' represents the contributions from the ith order approximation. 
The other independent variables are divided into i — I groups and arc arranged sequentially corresponding to the 
order of g-product decomposition. That is, the first variable group represents the first decomposition, the second 
variable group represents the second decomposition, and so on. Every variable group is a bit-string made of three 
possible element c, n, k and its length is equal to the number of the related order of g-product decomposition, that is, 
for the jth decompositions in the ith order approximation its length is i — j. In each variable group, c corresponds 
to a (5 function, n corresponds to a jy function and k corresponds to 1 (non-decomposition). Their sequence in the 
bit-string corresponds to the sequence of contraction and/or anti-contraction index string. From the above analysis 
and statement, the index string of the jth decompositions in the i order approximation is: 



n (7fc,7fc+i+j)- 



(55) 



fe=i 



For example, for A5, the first variable group is cccn, which refers to the first decomposition in five order approxi- 
mation and the terms to include the factor (57173(572 74*573 75 '774 76 ™ the definition of A^^cccn). Similarly, cncc means 
to insert the factor (57173777274^737557478 into the definition of A^lcncc). When there are nontrivial second contrac- 
tions, for instance, two variable groups {ccnn,kkc) represent that the definition of A^{cann,kkc) has the factor 
((57173 (57274 ?77375?77476) (5,^376- Since there are fully trivial contraction (the bit-string is made of only k), we omit their 
related variable group for simplicity. 

Furthermore, we pack up all the contraction- and non-contraction terms in the following way so that we can obtain 
conveniently the improved forms of perturbed solution of dynamics absorbing the partial contributions from the high 
order even all order approximations. We first decompose AJ^ , which is a summation of all above terms, into the 
three parts according to e~'^T«*, (-it)e~'^T«* and (— it)^e~'^Ti*/2: 



Ar' =Ar'{e) + A2'^'{te)+Ar'{t%). 



(56) 



Secondly, we decompose its every term into three parts according to e" 



lE-yt g lE^^t 



(E7, e-'^-i*) and e 



477 
477' 
477' 



(c) = A^'' (c-'^-'*) + A2^ (e-'^-i*) + AJ'' (e-'^v*), 
(te) = A^^'(te-'^-'*)+A^^'(te-'-^-i*)+A^^'(te-'^vt), 
it's) 



Al^ {fe-'^^*)+Al^ (i'e-'^-i*) + (i^e-'^v*). 



(57) 
(58) 
(59) 
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Finally, we again decompose every term in the above equations into the diagonal and off-diagonal parts about 7 and 



7' 



N), 



According to the above way, it is easy to obtain 



(60) 
(61) 
(62) 



where Ej. takes Ej,Ej-^ and Ej 



7ii72 

^r'(e-'^^*;N) = - Y^e 

71 



1 



(-E'7 ~ -^71) (-^7 ~ -^72)^ (-^7 ~ ^llf (-^7 ~ -^72) 



(i?'Y — -£^71 ) (-^7 — -^7' ) (^7 — -^71 ) (-^7 — -£/--y' ) 



^771^7172^727^^^,^ (63) 



„77l „7l7„77 

yi 5i yi 



7ii72 



g77 g7 7 g7 7 rjy^^ri^y' 

{E^ - Ej^ ) {E^ - E^2 ) {E^ - Eji ) ' 



(-£^7 ~ -^71) (-^71 ~ -^72) 



7li72 



: !Jl "tfl' '"tfi" ^ 'hn"h7' 

(£'7 ~ -^71) (-^71 ~ -^72) (-^71 ~ -^7') ' 



(64) 

(65) 
(66) 



Ar'(e-'^72t;D) = - ^ 



„77l „7l72 „727A 

i£^,t i/i i/i yi "77' 



e ^2 



71.72 



(E^ - Ej^) {E^^ - Ej^) 



7i>72 



77l„7l72 727 „ 
* t/l ffl ffl '/772'/77 

(-E'7 " ^72) (-E'71 ~ ^72) {^12 - -^7') 



^77'(e-i£yt.N) = ^e-'^v 

71 



+ 



(i/^ — £/-).') {Ej^ — ^7') (^7 ~ -^7') (^71 ~ -^7') 

ffl ^ g/ ^ ffl ' %i7'^77' 

(£^ - E^i) {E~f^ - Ey) {E^^ — Ey) 



„7 71 „7i7 „77 

9i 91 91 



(67) 
(68) 

(69) 
(70) 



In the end of this subsection, we would like to point out that the main purpose introducing the (7-product decom- 
position and calculating the contractions and anti-contractions of ^-product is to eliminate the apparent divergences 
and find out all the limitations from the contributions of g-product contraction terms. This is important to express 
the results with the physical significance. 



III. IMPROVED FORMS OF PERTURBED SOLUTION OF DYNAMICS 



In fact, the final aim using the ^-product decomposition and then calculating the limitation of the contraction terms 
is to absorb the partial contributions from the high order approximations of off-diagonal elements of the perturbing 

Hamiltonian matrix into the lower order approximations in our improved scheme of perturbation theory. In this 
section, making use of the skills and methods stated in previous section, we can obtain the zeroth, first, second and 
third order improved forms of perturbed solutions with the above features. 

In mathematics, the process to obtain the improved forms of perturbed solutions is a kind of technology to deal 
with an infinite series, that is, according to some principles and the general term form to rearrange those terms with 
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the same features together forming a group, then sum all of the terms in such a particular group that they become 
a compact function at a given precision, finally this infinite series is transformed into a new series form that directly 
relates to the studied problem. More concretely speaking, since we concern the system evolution with time t, we take 
those terms with (— iyit)e~'^'*, (— iyit)^e~'^'*/2! and (— i?/it)'^e~'^'*/3'j • • • with the same factor function / together 
forming a group, then sum them to obtain an exponential function /exp [— i (xj + y^) t\. The physical reason to do 
this is that such an exponential function represents the system evolution in theory and it has the obvious physical 
significance in the calculation of transition probability and perturbed energy. Through rearranging and summing, 
those terms with factors t°e~'^^>', (a = 1, 2, ■ • ■ ) in the higher order approximations are absorbed into the improved 
lower approximations, we thus can advance the precision, particularly, when the evolution time t is long enough. We 
can call it "dynamical rearrangement and summation" method. 



A. Improved form of the zeroth order perturbed solution of dynamics 



Let us start with the zeroth order perturbed solution of dynamics. In the usual perturbation theory, it is well-known 



(71) 



77' 



where a^i = (<&''' |\E'(0)). Now, we would like to improve it so that it can absorb the partial contributions from higher 
order approximations. Actually, we can find that ^2(c) and A^^nn^c) have the terms proportional to (—it) 



(-ii)e-'^-'' 
(-it)e-'^^* 



E 

L 71 

E 



1^771 |2 



^77' ' 



771 „7172 „727 

i/i 91 i/i 



Introduce the notation 



- E 

71 

G(3) . 



\9J 



1 



111 



71.72 



{Ey - EjJ{Ej - E^^) 



„77l „7l72^727 

91 9i 9i ■ 



(72) 
(73) 

(74) 
(75) 



It is clear that G^^^ has the energy dimension, and we will see that it can be called the ath revised energy. Let 
us add the terms ^ and the related terms in A4(te-'^^*, D), ^4(^26"*^^*^ D), ^5(^6"'^^*, D), A5(i^e~'^^*, D), 

A6(t^e~'^^*, D) and AQ{t^e~^^''*-) given in Appendix [Al together . that is. 



+ {-it) (g'^^) +0^^^) +0^^^^ +0^^^^ 
21 V 7 ^ ^7 



{-it) 



'^77' 1 



(76) 



Although we have not finished more calculations, from the mathematical symmetry and physical concept, we can 
think 



= e 



1 



(-it) (G(2)+G(3)+G(4)+Gf) 
+^ (g(2) + G'^^ + G(4) + Gf )' + . . 



'^77' ? 



(77) 



New terms added to the above equation ought to, we think, appear at A7(t), A%{t), Ai^{t) and Aio(t), or come from 
the point of view introducing the higher approximations. So we have 



and then obtain the improved form of the zeroth order perturbed solution of dynamics 



(78) 
(79) 



77' 
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(2) . (3) 

It is clear that G\ is real. In fact, G'y is also real. In order to prove it, we exchange the dummy indexes 71 and 
72 and take the complex conjugate of Gj , that is 



or = E 



igmgrngD 



E 

71.72 



(Ej - Ej^ ) {E^ - Ey^ ) 



„77l „7l72„72 7 

:9i ill 9i 



(80) 

where we have used the relations {jJi^^^'j ~ 9i^^^ for any /3i and (32 since Hi is Hermit. Similar analyses can be 
applied to G\ and G\ . These mean that e ^ -' "i ' is still an oscillatory factor. 



B. Improved form of the first order perturbed solution of dynamics 



Furthermore, in order to absorb the partial contributions from the approximation higher than zeroth order, we need 
to consider the contributions from off-diagonal elements in the higher order approximations. 
Well-known usual first order perturbing part of dynamics is 



|vi/(i)(o) = E 



-iE-,t 



Ey E-y' 



E-y — EyI 



-lE^t 



E^^ — E^' 



77 



(81) 



It must be emphasized that Hi is taken as only with the off-diagonal part for simplicity. That is, we have used the 
Hamiltonian redivision skill if the perturbing Hamiltonian matrix has the diagonal elements. 

Thus, from A^{te-'^-'\n) and A4(ie-'^^*, N), Ai{t^e-'^-'\V)), ^5(^26"'^^*, N), 746(^^6-'^^*, N), which are defined 
and calculated in the Appendix |^ it follows that 



AV it) 



[Ey - Ey ) 



(Ey - Ey) 

Therefore, by rewriting 



+ i-it) (g(2) + G(3) + G(4)) + ^ (Gf 

1 + , ^(3) , M^A , i-'ty 



' G\V + G^^ + G^f 



2! 



2, 2G(,2)G^3)^.. 



,-!(£-, +G<2)+G<3)+G<*))t -i(i5y+G<^' + G<^,' + G<1>)t' 



77 

gr 



77 

5i 



77 

5i 



(82) 



(83) 



we obtain the improved form of the first order perturbed solution of dynamics 

|*(i)(t))j = E^n'Wvl'i'^)- 

7,7' 



(84) 



C. Improved second order- and third order perturbed solution 

Likewise, it is not difficult to obtain 



= E - 



71 



,-i(iJ-,+G<^)+G<^))t _ ^i(iJ,,+G(^^)+G(^^))t 



{Ey - Ey^ 



9i gi O77' 



-!(£;-, +G!,^)+G<^))t 



A{E,,+G[f^+G[f^)t 



+ 



{Ej — -E-Yi) (-£^7 — ^7') (-^7 — ^71) (^71 ~ -^7') 

-i(^;,,+G<=,'4-G<^,')t 

_r aV^ gT^' > , 

{Ey Ery' ) {E'y-^ E'yf ) 



(85) 
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71,72 



-i(£;-,+G(^))t 



(^7 ^ -^71 ) (-^7 ~ -^72 )^ (^7 ^ ^'llf (-^7 ~ ^72 ) 



^7 71 / V 7 



^7 -^71/ V-^7 -^72 > 



(-E7 ~ ^llf (-^71 ~ -^72 ) (-^7 ~ -^72 f (^71 ~ -^72 ) 



E 

71 

E 



7 -^71/ v-^71 ^12) 

i{E-,+G'-^'')t 



-'1 ^11 1 V^7i -^72; 

-i(£;^+G<,^')t 



(-E-y — -E'l'i) (-E-y ^ (-^7 ~ -^71)^ (-^7 ~ -^7 



^771^7172^727^^^, 



„77l „7l7„77 

i/i yi i/i 



■^7 -"71 / V 7 7 
-i(£^+G<,^>)t 



'7772 



i(^^i+G?i')* 



^771 7' 



(-£^7 ^ -^71) (-^7 ~ -^72) (-^7 ~ -^7') (-^7 ~ -^71) (-^71 ~ -^72) (-^71 ~ -^7') 
^772 



g V '2 ' 72/7 



^771^7172^727 ^^^, 



(-^7 ~ ^72) (-E'71 - £^72) (-E'72 - £^7') 

Therefore, the improved forms of the second- and third order perturbed solutions are, respectively, 

7,7' 



(86) 

(87) 

(88) 



D. Summciry 

Obviously, our improved form of perturbed solution of dynamics up to the third order approximation is 

i=0 



(89) 



However, this solution absorbs the contributions from the whole Kf (to), A^^ (<^e) parts up to the fifth order 

approximation and the whole Ap (f^e), AJ^ (f^e) parts in the sixth order approximation. After considering the 
contractions and anti-contractions, we get the result corresponds to the replacement 



-iE^.t 



in the A^'^ (e) part, where 



(90) 



(91) 



i = 0,1,2, - ■■ , and 70 = 7- Here, we have absorbed the possible contributions from the diagonal elements of the 
perturbing Hamiltonian matrix. Although the upper bound of summation index a is different from the approximation 
order in the finished calculations, wc can conjecture that it may be taken to at least 5 based on the consideration 
from the physical concept and mathematical symmetry. For a > 5, their forms should be similar. From our point of 
view, such form is so delicate that its form happens impossibly by accident. Perhaps, there is a fundamental formula 
within it. Nevertheless, we have no idea of how to prove it strictly and generally at present. 

Actually, as soon as wc carry out further calculations, wc can absorb the contributions from higher order approxi- 
mations. Moreover, these calculations arc not difficult and are programmable because we only need to calculate the 
limitation and summation. Therefore, the advantages of our solution have been made clear in our improved forms 
of perturbed solution of dynamics. In other words, they offer clear evidences to show our improved scheme is better 
than the existing method in the precision and efficiency. In the following several sections, we will clearly demonstrate 
these problems. 



IV. IMPROVED TRANSITION PROBABILITY AND REVISED FERMI'S GOLDEN RULE 



One of the interesting applications of our perturbed solution is the calculation of transition probability in general 
quantum systems independent of time. It ameliorates the well-known conclusion because our solution absorbs the 
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contributions from the high order approximations of the perturbing Hamiltonian. Moreover, in terms of our improved 
forms of perturbed solution, it is easy to obtain the high order transition probabihty. In addition, for the case of 
sudden perturbation, our scheme is also suitable. 

Let us start with the following perturbing expansion of state evolution with time t, 



7 n—0 7 



(92) 



When we take the initial state as |^^), from our improved first order perturbed solution, we immediately obtain 



.(1) 

"7,1 



E~, — Ea 



5i I 



(93) 



where 



(94) 



Here, we use the subscript "I" for distinguishing it from the usual result. Omitting a unimportant phase factor e '^t*, 
we can rewrite it as 



7/3 

9i 



''■^ - Ei3 



(95) 



where uj^p — E^ — Efj. Obviously it is different from the well known conclusion 



(96) 



where w^/j = Ey ~ Ep. Therefore, our result contains the partial contributions from the high order approximations. 
Considering the transition probability from \^^) to \(^'^) after time T, we have 



7/3 



1 _ gl"7|9J 



7/3 



sir? (w7/3r/2) 

(^7/3/2)^ 



(97) 



In terms of the relation 



sin^ X — SIT? y = — [cos(2y) — cos(2a;)] 



we have the revision part of transition probability 



7/3 

9{ 



) (ujjfjT) — cos (w-y/jT) 



(99) 



If plotting 



{U^p/2y V^7/3/ (W7/3/2) 



2 „• 2 



sin^ (^7/3^/2) _ f^-ypV sin^ {^'ypT/2) 



(100) 



we can see that it has a well-defined peak centered at uj.yp — 0. Just as what has been done in the usual case, we can 
extend the integral range as —00 —^ 00. Thus, the revised Fermi's golden rule 



w — wf+ Aw, 



(101) 



where the usual Fermi's golden rule is 



WF = 2TTp{Ep) 



7/3 

9i 



(102) 
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in which w means the transition velocity, p{E-y) is the density of final state and we have used the integral formula 



°° sin^x 



while the revision part is 



Aw ^2 dE^p (E^) 



7/3 

9{ 



cos (ujjfsT) — cos {uj-yfjT) 



(103) 



(104) 



It is clear that is a function of E-^, and then a function of uj-fp. For simplicity, we only consider uj^p to its second 
order approximation, that is 



\9T 



71 

Again based on dE^ = dw-y/j, we have 

/>oo 

Aw = 2 I duj^pp {uj^p + Ep) 



^7/3 — 'i'71/3 tt^/371 



0{Hl). 



(105) 



7/9 

9T 



[ijj^pT] — cos [uj-yp (w) T] 



(106) 



It seems to not to be easy to deduce the general form of this integral. In order to simplify it, we can use the fact that 
U)^I3 — Lu^i3 is a smaller quantity since 



G 



(107) 



For example, we can approximatively take 



cos (ujjpT) — cos {oj^pT) « T {uj-yp — uj-yp) sin {uj^pT — uj^pT) , 



(108) 



then calculate the integral. We will study it in our other manuscript (in preparing). 

Obviously, the revision comes from the contributions of high order approximations. The physical effect resulted 
from our solution, whether is important or unimportant, should be investigated in some concrete quantum systems. 
Recently, we reconsider the transition probability and perturbed energy for a Hydrogen atom in a constant magnetic 
field [Sj]. We find the results obtained by using our improved scheme are indeed more satisfying in the calculation 
precision and efficiency. We will discuss more examples in our future manuscripts (in preparing). 

It is clear that the relevant results can be obtained from the usual results via replacing Lj^p in the exponential 
power by using ui-yp. Thus, one thing is true — with the time t evolving, e^'^^T^^"*/^^ term in the improved transition 
probability can be very different from e^'^"^''*/^' in the traditional one, which might lead to totally different results. 
To save the space, we do not intend to discuss more here. 

In fact, there is no any difficulty to obtain the second- and three order transition probability in terms of our 
improved forms of perturbed solution in the previous section. More higher order transition probability can be given 
effectively and accurately by our scheme. 



IMPROVED FORMS OF PERTURBED ENERGY AND PERTURBED STATE 



Now we study how to calculate the improved forms of perturbed energy and perturbed state. For simplicity, we 
only study them concerning the improved second order approximation. Based on the experience from the skill one 
in Sec. IVIl we can, in fact, set a new E and then use the technology in the usual perturbative theory. That is, we 
denote 



E^^ — E^^ 



(109) 
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7,7' ^. 



E 

71 



71 



2 i/l Hi "77 



(-E^ — -E'71 ) (-£^7 — -E^7') (-^7 ^ ^71) (-^71 ~ -^7') 



J7 ij7j y (^jj7 



(110) 



Note that £'7. can contain the diagonal element /i7' of the original Hi, and we do not obviously write /i]'' and take 
new Hi matrix as ofF-diagonal in the Ho representation. 
Because that 



(111) 



we have 



7' 



E^ Ery' 77' 

E^ — Ej' 



E ^7 -^71 „77l „717 X , 
— —r29i 9i Oil' + 2^ 
71 (Ef-Ey-,) 



E^ 



(£'7 — £71 ) (-£7 — £7' 



71 



£. 



■7' 



771 717 I 

W 5l '777' } 



{Ej — £71) (£71 ~ -^7') (^7 ~ -^7') (^71 — ^7') 
In the usual perturbation theory, Hi is taken as a perturbing part with the form 

Hi = Xv, 



(112) 



(113) 



where A is a real number that is called the perturbing parameter. It must be emphasized that E^. can be taken as 
explicitly independent perturbing parameter A, because we introduce A as a formal multiplier after redefinition. In 
other words, E^. has absorbed those terms adding to it and formed a new quantity. This way has been seen in our 
Hamiltonian redivision skill. Without loss of generality, we further take Hi only with the off-diagonal form, that is 



^^7172 =^7172 ^ ;^y7l72^ 



(114) 



Then, we expand both the desired expansion coefficients a-y and the energy eigenvalues Et in a power series of 
perturbation parameter A: 



Et — E-^'-^T,!' 0^7 ~ E "^'"-v'r 



(115) 



1=0 



1=0 



1. Improved 0th approximation 

If we set A = 0, eq. (fTT2)) yields 

E^\a^l^E,a^l (116) 
where 7 runs over all levels. Actually, let us focus on the level 7 = /3, then 

4°i = ^P- (117) 

When the initial state is taken as | 

= S,0. (118) 

Obviously, the improved form of perturbed energy is different from the results in the usual perturbative theory because 
it absorbs the contributions from the higher order approximations. However, the so-call improved form of perturbed 
state is the same as the usual result. 
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2. Improved 1st approximation 

Again from eq. pl2l) it follows that 



7 ' ' 



When 7 = /3, it is easy to obtain 



If 7 ^ /3, then 



i;^^] = 0. (120) 



It is clear that the first order results are the same as the one in the usual perturbative theory. 

3. Improved 2nd approximation 

Likewise, the following equation 

+ i;(i)a(i) + ^;(o)a(2) _ (2) , ^7 - -^7' 77' (i) _ -^l^ll^y-yii^nT 5 ,^(0) 

7' T '''' 71,7' 1^7 "^71/ 



71,7' L 



{E^ — -£'71) (^7 ^ -^7') (-^7 — -^71) (^71 ^ ^7') (-^7 ^ -^7') (-^71 — ^7') 
is obtained and it yields 



^77i„7i7^^^,aW^. (122) 



E^^^^ = 0, (123) 

if we take 'y = f3. When 7 7^ /3, we have 

It is consistent with the off-diagonal part of usual result. In fact, since we have taken Hj'^ to be off-diagonal, it 
does not have a diagonal part. However, we think its form is more appropriate. In addition, we do not consider the 
revision part introduced by normalization. While Ej,.-^ = is a new result. 

4. Summary 

Now we can see, up to the improved second order approximation: 

^Ep^E0 + Gf + Gf . (125) 

Compared with the usual one, they are consistent at the former two orders. It is not strange since the physical law is 
the same. However, our improved form of perturbed energy contains a third order term. In other words, our solution 
might be effective in order to obtain the contribution from high order approximations. The possible physical reason 
is that a redefined form of the solution is obtained. 

In special, when we allow the H'J^ to have the diagonal elements, the improved second order approximation becomes 

Et,p -Ep + hl + Gf + Gf . (126) 
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Likewise, if we redefine 

E^, = E^, + hf + + G^3) ^ Qi4) _ (127) 
Thus, only considering the first order approximation, we can obtain 

Et,0 ^E0 + h^,+ of + of + of . (128) 
In fact, the reason is our conjecture in the previous section. The correct form of redefined E^^ should be 

Et.j3 ^E(, + h^( + Gf + Gf + Gf + + • • • . (129) 

This implies that our improved scheme absorbs the partial even whole significant contributions from the high order 
approximations. In addition, based on the fact that the improved second approximation is actually zero, it is possible 
that this implies our solution will fade down more rapidly than the solution in the usual pcrturbativc theory. 

Actually, the main advantage of our solution is in dynamical development. The contributions from the high order 
approximations play more important roles in the relevant physical problems such as the entanglement dynamics and 
decoherence process. For the improved perturbed energy, its high order part has obvious physical meaning. But, 
for the improved form of perturbed state, we find them to be the same as the existed perturbation theory up to the 
second approximation. 



VI. EXAMPLE AND APPLICATION 



In order to concretely illustrate that our exact solution and the improved scheme of perturbation theory are indeed 
more effective and more accurate, let us study a simple example: two-state system, which appears in the most of 
quantum mechanics textbooks. Its Hamiltonian can be written as 

"'{It)- w 

where we have used the the basis formed by the unperturbed energy eigenvectors, that is 

i*')=(i)' i*')=(?)- ^'''^ 

In other words: 

Ho\^^)=E^\^^), (7 = 1,2) (132) 



where 



Thus, this means the perturbing Hamiltonian is taken as 

^. = (y°/fl- (134) 



This two state system has the following eigen equation 



7 

It is easy to obtain its solution: corresponding eigenvectors and eigenvalues 



= E'£\^^). (135) 



= 1 ( + ^li \ n36) 
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(138) 



where \V\ — IV12I = IV21I, UJ21 = E2 — Ei, LU21 — E2 — E^ = \/'i\V\'^ + uj^-^, and we have set E2 > Ei without loss of 
generahty. 

Obviously the transition probability from state 1 to state 2 is 



P^(1^2) = |($2|e-'«*|$i)|' = 



71.72 = 1 



_ sin^K,t/2) 



In the usual perturbation theory, up to the second order approximation, the well-known perturbed energies are 

12 |y|2 



E[ = E, 



' -, E^ — El 



UJ21 LU21 

While, under the first order approximation, the transition probability from state 1 to state 2 is 



(140) 



(c^2l/2)2 • 

Using our improved scheme, only to the first order approximation, we get the corresponding perturbed energies 



141) 



iil — £/! 1 , ii2 — J^l 



UJ21 UJ21 



UJ21 UJ21 



where we have used the facts that 



/^(2) _ I^P _ ^(2) ^(3) _ ^(3) _ „ ^(4) _ \V\^ _ ^(4) 
UJ21 id2i 



(142) 



(143) 



Obviously, under the first order approximation, our scheme yields the transition probability from state 1 to state 2 as 

jsin^ (w2it/2) 



Pi(1^2) = |y|^ . 

(W2l/2)^ 

where tD2i = E2 — Ei. Therefore we can say our scheme is more effective. Moreover, we notice that 



Ei = El 



\V\' , \V\ 



\v\' 



+ ^ + o{\vn ^Ei + o{\vn = E[ + ^ + oi\vn 



Ei 



El 



and 



UJ21 UJ21 

\v^jvl_ 

UJ2I ^21 



^21 



oi\vn. 



P^(l - 2) = \V\ 



2 sin^_(tJ2ii/2) 

(cc'2l/2)2 



sin {u}2it) sin^ (0^21^/2) 
2(l^2i/2)2 " (a;2i/2)3 



{UJ21 - W21) + 0[{UJ21 ~ UJ2iy 



, ,nSin^ (w2ii/2) 
Pi(l ^ 2) - \V\^ V 21 / ^ 



{UJ21 - W21) + ©[(^21 - UJ2l)'^] 



= P(l ^ 2) + \V\ 



(a;2i/2)3 
sin(ti)2i^) sin^ (aj2ii/2) 



_2(c^2i/2)2 (w2i/2)3 
Therefore, we can say that our scheme is more accurate. 



(tJ21 - tJ2l) + ©[(^^21 - t^2l) 



(144) 

(145) 
(146) 

(147) 
(148) 
(149) 



VII. DISCUSSION AND CONCLUSION 



In this paper, our improved scheme of perturbation theory is proposed based on our exact solution in general 
quantum systems Q. Because our exact solution has a general term that is a c- number function and proportional 
to power of the perturbing Hamiltonian, this provides the probability considering the partial contributions from the 
high order even all order approximations. While our dynamical rearrangement and summation method helps us to 
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realize this probability. Just as the contributions from the high order even all order approximations are absorbed to 
the lower order approximations, our scheme becomes an improved one. 

It must be emphasized that our improved scheme of perturbation theory is proposed largely dependent on the facts 
that we find and develop a series of skills and methods. From that the Hamiltonian redivision skill overcomes the 
flaw in the usual perturbation theory, improves the calculation precision, extends the applicable range and removes 
the possible degeneracies to that the perturbing Hamiltonian matrix product decomposition method separates the 
contraction terms and anti-contraction terms, eliminates the apparent divergences in the power series of the perturbing 
Hamiltonian and provides the groundwork of "dynamical rearrangement and summation" , we have seen these ideas, 
skills and methods to be very useful. 

Actually, the start point that delays to introduce the perturbing parameter as possible plays an important even key 
role in our improved scheme of perturbation theory. It enlightens us to seek for the above skills and methods. 

From our exact solution transferring to our improved scheme of perturbation theory we does not directly use the 
cut-off approximation, but first deals with the power series of perturbation so that the contributions from some high 
order even all order approximations can be absorbed into the lower orders than the cut-off order as possible. Hence, our 
improved scheme of perturbation theory is physically reasonable, mathematically clear and methodologically skillful. 
This provides the guarantee achieving high efficiency and high precision. Through finding the improved forms of 
perturbed solutions of dynamics, we generally demonstrate this conclusion. Furthermore, we prove the correctness 
of this conclusion via calculating the improved form of transition probability, perturbed energy and perturbed state. 
Specially, we obtain the revised Fermi's golden rule. Moreover, we illustrate the advantages of our improved scheme 
in an easy understanding example of two-state system. All of this implies the physical reasons and evidences why our 
improved scheme of perturbation theory is actually calculable, operationally efficient, conclusively more accurate. 

From the features of our improved scheme, we believe that it will have interesting applications in the calculation of 
entanglement dynamics and decoherence process as well as the other physical quantities dependent on the expanding 
coefficients. 

In fact, a given lower order approximation of improved form of the perturbed solution absorbing the partial con- 
tributions from the higher order even all order approximations is obtained by our dynamical rearrangement and 
summation method, just like "Fynmann figures" summation" that has been done in the quantum field theory. It 
is emphasized that these contributions have to be significant in physics. Considering time evolution form is our 
physical ideas and absorbing the high order approximations with the factors i°e~''^"'i*, (a = 1, 2, ■ • ■ ) to the improved 
lower order approximations definitely can advance the precision. Therefore, using our dynamical rearrangement and 
summation method is appropriate and reasonable in our view. 

For a concrete example, except for some technological and calculational works, it needs the extensive physical 
background knowledge to account for the significance of related results. That is, since the differences of the related 
conclusions between our improved scheme and the usual perturbation theory are in high order approximation parts, 
we have to study the revisions (differences) to find out whether they are important or unimportant to the studied 
problems. In addition, our conjecture about the perturbed energy is based on physical symmetry and mathematical 
consideration, it is still open at the strict sense. As to the degenerate cases including specially, vanishing the off- 
diagonal element of the perturbing Hamiltonian matrix between any two degenerate levels, we have discussed how 
to deal with them, except for the very complicated cases that the degeneracy can not be completely removed by 
the diagonalization of the degenerate subspaces trick and the Hamiltonian redivision skill as well as the off-diagonal 
element of the perturbing Hamiltonian matrix between any two degenerate levels are not vanishing when the remained 
degeneracies are allowed. 

It must be emphasized that the study on the time evolution operator plays a central role in quantum dynamics 
and perturbation theory. Because of the universal significance of our general and explicit expression of the time 
evolution operator, we wish that it will have more applications in quantum theory. Besides the above studies through 
the perturbative method, it is more interesting to apply our exact solution to the formalization study of quantum 
dynamics in order to further and more powerfully show the advantages of our exact solution. 

In summary, our results can be thought of as theoretical developments of perturbation theory, and they arc helpful 
for understanding the theory of quantum mechanics and providing some powerful tools for the calculation of physical 
quantities in general quantum systems. Together with our exact solution 3] and open system dynamics [6], they 
can finally form the foundation of theoretical formulism of quantum mechanics in general quantum systems. Further 
study on quantum mechanics of general quantum systems is on progressing. 
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APPENDIX A: THE CALCULATIONS OF THE HIGH ORDER TERMS 



Since we have taken the Hi only with the off-diagonal part, it is enough to calculate the contributions from them, fn 
Sec. I VII the contributions from the first, second and third order approximations have been given. In this appendix, we 
would like to find the contributions from the fourth to the sixth order approximations. The calculational technologies 
used by us are mainly to the limit process, dummy index changing and summation, as well as the replacement 
Qi^^Viiij = 9i'^'^ since g'^^^ has been off-diagonal. These calculations are not difficult, but are a little lengthy. 



1. 1 = 4: case 



For the fourth order approximation, its contributions from the first decompositions consists of eight terms: 

AJ*' = A}''' (ccc) + A}^' (ccn) -t- AJ^' (cnc) + AJ'' (ncc) 

+Aj''''{cnn) + A'f'incn) + A'f'innc) + Aj'''{nnn). (AI) 

Its the former four terms have no the nontrivial second contractions, and its the fifth and seven terms have one 
nontrivial second contraction as follows 



AJ'' (cnn) — AJ'' {cnn, kc) + AJ'' {cnn;kn), 
A'P {ncn) ~ AJ'^ {ncn,c) + AJ'^ {ncn,n), 
A'P (nnc) — AJ'^ [nnc, ck) + AJ'^ {nnc,nk). 



(A2) 
(A3) 
(A4) 



In addition, the last term in eq. (|AI|) has two nontrivial second contractions, and its fourth term has also the third 
contraction. Hence 



AJ^ (nnn) — AJ^ (nnn, cc) + AJ'^ {nnn,cn)A']^ [nnn^nc) + A'^^ [nnn,nn), 
AJ^ {nnn,nn) = AJ^ (nnn,nn,c) + AJ'^ {nnn,nn,n). 



(A5) 
(A6) 



All together, we have the fifteen terms that are the contributions from whole contractions and anti-contractions of 
the fourth order approximation. 

First, let us calculate the former four terms only with the first contractions and anti-contractions, that is, with 
more than two S functions (or less than two rj functions) 



AT' (ccc) = 



71 



.4=1 



3e 



3e-'^^i* 



4 

n 



9i 



2e 



n '^7fc7fc + 2 1 '^717*^757' 
\fe=l / 



-(-it) 



^7 -^11 

{Ey - E^^ ) 



^7 



2 {E^ - Ej^)'^ 



(A7) 



AJ-" (c 



:75 Li— I 



E 
E 



E(-ir^ 



■J~J „7j7j + l 



9i 



'^7173 '^7274^7375 '^7 17 "^757' 
2g-iB^t 



{E^ — £"71) (-£'7 ~ E-y') [E-y — -E'71)'^ {Ef — E^i) {E^ — ) (-£'71 — -£7') 



J7 JJ7J j yi^^ J^y' 



(E-y ~ Ej, f {E^, - Ey ) {E-y - Ey f {E. 



V7 -1J71; y-^ti -L^i 



71 -^7' 



+ (-it) 



Q—iEjt 



(i?7 — -£71 ) (-£7 — ) 



(£'7 — i?7j ) (-£71 — E^' 



(A8) 
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AT {cnc) = 

E 



7ir-- >75 Li=l 



n 



„7j73 + l 

yi 



'^71 73 %2 74 "^73 75 '^71 7 '^75 7' 



7l>72 



+ 



Q—iE^t 



{E-y - Ej^) {E^ - E^^)^ [E^ - E^^Y [E^ - E^^Y {E^ - Ej^Y {E^ — Ej^) 



e "2 



(-E'7 - -^71 )^ (-^71 ~ -^72 ) (-^7 ~ -^72 )^ (-^71 ~ 



i-it) 



-i-it) 



7 72 ^ V-^71 72 } 



{E^ -E^,Y {E^ - ^72 ) 2 (i?^ - i?^, ) {E^ - E^, ) 
11' ■ 



(A9) 



AT\ncc) = 

IT' 

= E 



7l)---)75 Li=l 



4 

n 



„73 7j + l 



^7173*^7274 '^7375 ^717^757' 



(i^/y — -^71 ) (-^7 — ) (-^7 — -^71 ) {E^i — E'^f ) (-^7 — -^71 ) {E^i — E^f ) 



(iz-y — E^') {Ej-^ — E-y'Y {E'y ~ E-y'Y {E~fi — E-i'Y {E-y — ErfY) (-^71 — E^'Y 



-i-it) 



^717 



grgr ^77' 



(AlO) 



Then, we calculate the three terms with the single first contraction, that is, with one 5 function. Because one 5 
function can not eliminate the whole apparent singularity, we also need to find out the nontrivial second contraction- 
and/or anti-contraction terms. 



A]Y {cnn, kc) = 

71, •• 

= E 



E(-ir^ 



di{E[j,4]) 



4 

n 



„73 7j + l 

9i 



^71 73 "^12 74 ^773 75 '^72 7' '^71 7 <57 



2g-iB.yt 



Q—iE^t 



+ 



Q—iE'y^ t 



{E'y — -^71) {-^7 — -^7')^ (-^7 — -^71)^ (-^7 — -^7')^ (-^7 — -^71)^ (-^71 — -^7')'^ 



2c 



-lE.t 



{E^ — -^7') (-^71 — E'yf) {E^ — E'yf) (-^71 — Ejf) 

2 



+ (-ii) 



^7 -^7'y V-^71 -^7 
-(-ii) 



g-iByt 



(iz-y — -£'7') (-£71 -^7') 



(£^7 — -^71) (-^7 — -^7') 



77 

5i 



^771^717 



(All) 



^42''' {cnn,kn) 



E 

7l>--- >75 

E 

7li72 



E(-ir^ 



MEYA] 



.9i 



•^71 73 ^72 74 ^73 75 ^72 7' ^71 7 ^' 



'757' 



(-E'7 ~ -^71) (-^7 ~ -^72) (-^7 ~ -^7')^ {Ej - E^Y) (-^7 ~ -^72)^ (-^7 ~ -^7') 



g — 1^7* 



g — 1^71 1 



(-£^7 ~ -^71)^ {E^ - E^Y {E-i - E^i) {E^ - E^Y^ {Ej^ - E^Y (-^71 ~" -^7') 



e~'^72* 



+ 



g-ii5yt 



(-£^7 ~ -^72)^ (-^71 ~ -^72) (-^72 ~ E^i) {Ej - EyY {Eji - Ey) {Ej^ - Ej') 



+(-«) 



Q—iE^t 



(-^7 ~ -^71) (-^7 -^72) {Ej - Ey) 



1 9i^^ I 9i"'^ ^71 72 ^71 7' ^77' 



(A12) 



24 



AP [ncn, c) 



7ir-- 175 Li=l 

E 

7li72 



n 



„73 7j + l 

91 



Q—iE^t 



'77l73'^7274'?7375'^7l7^757''^77 
2g-iE^t 



Q—iE^^ t 



{E^ - E^^f [E^ - Ej^f [E^ - E^^ )^ {Ej - E^^ ) {E^ - E^if (-^71 ~ -^72 )^ 



{Ej - Ej^ f {E^^ - Ej^ ) {E^ - E^^ f {E^^ - E^^ f {E^ - E^if {Ej - E^^ ) 



(A13) 



Ap {ncn,n) 



E 

It'" 

= E 



E(-ir^ 

7l;--- .75 Lj=l 



4 

n 



„73 7j + l 

9i 



Vji 73 ^72 74 ^7375 '^717 "^75 7 ' '?77' 



(-E'7 ~ -^71)^ (^7 ~ -^72) (-£^7 ~ -^7') (-^7 ~ -^71) (-^71 ~ -^72) (-£'71 ~ -^7')^ 



+ 



g— i-E^i * 



g— i-E^i* 



{Ej - E^^) {Ej^ - E-f^f {E^^ - Ey) {E^ - E^if (^71 ~ -^72) (-^71 ~ ^7') 

g-i£^2* 



+ ■ 



^7 ^11 J -^72 J V-^71 -^7 

g-iByt 



(■£^7 - -£^72) (■£'71 ~ -£'72)^ (-£'72 ~ -£'7') (-£'7 ~ -£'7') (-£'71 ~ -£'7')^ (-£'72 ~ -£'7') 



-{-it) 



{Ej - E^^) {Ej^ - E^^) {Ej^ - Ej>) 



|^7i72| gj^^gj^'^ '7772 ^727'%7'- 



(A14) 



{nnc,ck) = 



7ij--- i76 Lj=l 



E(-ir^ 



-iE-, t 



4 

n 



„73 7j-I-i 

9i 



'77173^7274^7375 ^71 74 ^717 ^757' 



= E 



— -^71) (-^7 — -^7')"^ (^7 — -^71) (^7 — -^7') (-^7 — -^71) (-^71 — -^7') 



71 L 

e-'^y* 

(S-y — -E-y') (-E71 — -£'7') (-£^7 — -£'7')"^ (-£'71 ~ -£'7') 



^7 -^71 J V-^7 -^7' 
2g-i£y* 



+ (-it) 



(-E-y — -^71) (-^7 ~ -^7') 



-(-it) 



Q-iE^t 



{E'y E^' ) (-^71 -E-j,' ) 



^771^717 



(A15) 



A]^ {nnc, nk) 



E 

7ir-- 175 

E 

71.72 

+ 



E(-ir' 



rfi(i^[7,4]) 



7j7j + i 



'771 73 Vj2l4 ^7375 '77174 '^717^767' 



g — 1-^71* 



(£y - E'yJ (iJ-y " Ej^) {E^ " -B^')^ (£7 " EjJ {E^^ - E^^) {E^^ - E^)^ 



g-iB^jt 



g-iByt 



(-£^7 ~ -£'72) (-£'71 ~ -£'72) (-£'72 - -£'7')^ (-£^7 ~ -£'7') (-£'71 ~ -£'7') (-£'72 - -£-7')^ 

g-iSy* 



■^71 -^7 ) V-^72 
g-iByt 



(-Ey — -£'7') (-£'71 — -£'7')'^ (-£'72 — -£^7') (-£'7 ~ -£'7')^ (-£^71 — -£^7') (-£'72 — -£'7') 

-{-it) 



g-iByt 



(-£^7 ~ -£'7') (-£'71 ~ -£'7') (-£'72 ~ -£'7' 



7 -^7 / V-^71 -^7 M-^72 -^7 
2 



^727 



^77 gl 7 7^^^2?77i 72^777' 



(A16) 



Finally, we calculate the {nnn) by considering the two second decompositions, that is, its former three terms 
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AJ'^ {nnn,cc) = 

71, ■• 

= E 



7l,--- i75 Li=l 



4 

n 



„737j + i 
i/1 



'niklk+2 1 ^7174 '^7275 '^717^' 



'757' 



2c 



Vfc=l 

-li"-, / 



(£/-y — E'yi) (-£'7 — -^7')*^ (-^7 — -^71)'^ (-^7 — -^7')^ (-^7 — -^71)^ (-^71 — -^7')'^ 



2e-'^7 



(i/-y -^7') (-^71 -^7') (^7 -^7')^ (-^71 -^7') 



+ (-ii) 



— -^71) (-^7 — -^7') 



-(-it) 



5 



„77 „7 71 „7i7„77 



A4''' {nnn,cn) 



E 

7ir-- >• 



4 

n 



„7j73 + l 
ill 



Y\ ^7fc7fc+2 1 ^7174 ^72 76 '^71 7 '^75 7' 



g— iS-,t 



(£'7 - -S71) (-^7 ~ -^72) (-^7 ~ -£^7')^ (-E'7 - -E71) (-E7 - -£72)^ (£'7 - -E'7') 



{E^ - -^71)^ (£'7 - -E72) (£7 - £7') (£7 - £71)^ (-E'71 - £72) (£71 - £7') 



■^71/ V-^71 ^12 J V-^71 7 



(£7 ^ £72)^ (£71 ~ £72) (£72 - £7') (£7 ~ £7')^ (£71 ~ £7') (£72 - £7') 



AJ^ {nnn,nc) 



7ir-- >75 



E 

1 )"■ 

E 



(£7 ~ £71) (£7 ~ £72) (£7 ~ £7') 

El-')-*" 



.(£[7,4]) 



m'j+i 



JJ[ Vlklk+2 1 '77l74^7275^7l7'^757' 



^—'lE^t 



(£7 ~ £71) (£7 ~ £72) (£7 ~ £7')^ (£7 ~ £71) (£71 ~ £72) (£71 ~ £7')^ 



+ 



(£7 - £72) (£71 ~ £72) (£72 ~ £7')^ (£7 ~ £7') (£71 ~ £7') (£72 ~ £7')^ 



e-iEyt 



(£7 ^ £7') (£71 ~ £7') (£72 - £7') (£7 ~ £7') (£71 ~ £7') (£72 - £7') 

-(-ii)- 



(£7 - £7') (-^71 ~ -^7') (-^72 - £7') 

while the fourth term has the third decomposition, that is 



^77 71^7172^727 rj^^^rjry^^. 



(A17) 



(A18) 



(A19) 



A]^ {nnn,nn,c) = 

7i>--- >7 

= E 



E(-ir^ 



i=l 



^i.(£[7,4]) 



J|^^7,7. + l 



fc7fc+2 (57i7<5757'577' 



Q—iE^t 



(£7 ~ £71) (£7 ~ £72) (£7 ~ £73)^ (£7 ~ £71) (£7 ~ £72)^ (£7 ~ £73) 



Q—iE^t 



+ 



Q—iEj^ t 



(£7 ~ £71)^ (£7 ~ -^72) (£7 ~ £73) (£7 ~ £71)^ (£71 ~ £72) (£71 " £7 



Q-iEy^t 



+ 



(£7 ~ £72)^ (£71 ~ £72) (£72 - £73) (£7 ~ £73)^ (£71 ~ -^73) (£72 - £73) 



^7 -^72/ V-^71 -^12 J V-^72 -^73 
Q—iEyt 

^ {E^ - E^^) {E^ - E^^) {E^ - £'^3) 



■'1 -^73/ V-^71 -^73; V-^72 -^73; 

^771 „7i72 _7273 „737» 



77i„7i72 7273„737„ „ X 
i/l i/l i/l ai '/772'/7i73"' 



77 



(A20) 



26 



nnn, nn, n] 



E 



,75 Li=l 

E 



n 9^'^' 



(i?-y - Ej^ ) (£;^ - ) {Ey - Ey^ ) {E^ ~ Ey ) 

a-iE^, t 



{Ey - E^J (Ey-^ - Ej^) {E-y^ - Ey^) {E^^ - Ey) 

-'iE~.^ t 



{Ey - E-y^) {Ej-^ - E^^) (£'^2 - E^^) [E^^ - £;^/) 

[E^ - £^73) (-£71 ~ -£73) (-E'72 - -5^73 ) (-£73 - £7') 
(i?^ - E^i) {E~f^ - Ey) {E^^ - Ey) {E~f,^ - Ey) 



^9l 9l 9l 9l '7772^773'?77''?7l73^7l7''?727'- 

(A21) 

Now, all 15 contractions and/or anti-contractions in the fourth order approximation have been calculated out. 

In order to absorb the contributions from the fourth order approximation to the improved forms of lower order 
perturbed solutions, we first decompose AJ^ , which is a summation of all above terms, into the three parts according 
to their factor forms in e"'^^;*, (— i<)e~'^^i* and {—it)^e~'^^i* /2, that is 



A^ =Al''ie) + Ar{te)+Arife). 



(A22) 

-i-E-,, t\ 



Secondly, we decompose its every term into three parts according to the factor forms in e '^T*,e (^^^ e 

and e"'^T'*, that is 



(A23) 
(A24) 
(A25) 

Finally, we again decompose every term in the above equations into the diagonal and off-diagonal parts about 7 and 
7', that is 



AJ-* (e) - AJ'' (e-'^-*) -t- (e-'^-i*) + AJ'' (e-^v*), 
Aj'^'ite) = A2^'(te-'^-'*)+A2^'(te-'^-'i*)+^r'(te-'^-'*), 



^77 (g-iis,.t) ^ yl]^ (e-'^-.*;D)+A2^ (e-'^-.*;N), 



(A26) 
(A27) 
(A28) 

where E-y. takes E^ , E^-^^ and Ejr . 

If we do not concern the improved forms of perturbed solutions equal to or higer than the fourth order one, we 
only need to write down the second and third terms in ea. (|A22p and calculate their diagonal and off-diagonal parts 
respectively. Based on the calculated results above, it is easy to obtain 



/1 77 
^4 



-iE^t, 



D) = (-ii)e 



-iE^t 



E 



-2 15/ 



771 1 



71 (Ey - Ej^) 

j2 ig7"i'igri'^7i72 

71,72 (-^7 ~ -^71) (£7 ~ £72) 

- E 

71,72,73 



E 



7 -£72 J 



71,72 (-^T ~ -^Tl) i^' 
71,72 (^'T ~ ^11) (-£7 ~ £72) 



^7772 



„77l„7l72 7273 737„ „ 
91 9 1 t/l '/772'/7l73 

(-E'7 ~ £71) (-^7 ~ ^72) (-£7 ~ -^73) 



(A29) 



Substituting the relation riji^pr^ = 1 
of so-called the ath revision energy G 



0/31/32: 



using the technology of index exchanging and introducing the definitions 



- E 



71 



\9V'\" 



(A30) 
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„77l „7l72 „7273 „737„ 
ffl t/l j/l t/l '/772 

{Ej - Ej^ ) - Ej2 ) (£'^ - Ey^ ) 



71.72,73 

we can simplify Eq. (|A29[) to the following concise form: 

(ie-'^-*;D) = -(~it)e-'^-* 
Similar calculation and simplification lead to 

yir' (te-'^-*;N) = (-it)e~'^-* 

where 



„77l „7l7„772 „727 
ffl ffl ffl t/l 

71,72 (-^■T ~ ■^'Tl) (^7 ~ -^72 



E 



E 



G 



(2) 



(■£^7 ~ -^71 ) 



L 71 ^ T T 



+ E 



(2) 



/^(3)„77' 

(-E-y — Ej/) (-^7 — -^71 ) (-^7 ^ -^7') 



„77l „7l72^727 

^(3) ^ ^ 9\ 9\ 9i 



71 ,72 



(-^7 ~ -^71 ) (^7 ~ ^72) 



(A31) 



(A32) 



(A33) 



(A34) 



For saving the space, the corresponding detail is omitted. In fact, it is not difficult, but it is necessary to be careful 
enough, specially in the cases of higher order approximations. 
In the same way, we can obtain: 



(-i^E 



2(2)Q-i_E^, t 
71 

{Ej - Ej-, 



77 ; 



G^'^e-'^-i* 



71 



(^7 ~ -^71 ) (^71 ~ -^7' ) 



771 7l7 

Si Si VfY 



(A35) 
(A36) 



A2-''(ie-'^y*;D) = 0, 

^2^' (te-'^-'*;N) = -(-ii)e-'-^vt 

+E 



E 



^(3) 

*^7' 77' 

— :9i 



71 



(2) 



g;i 



71 



{Ej — Ej') {Ej-^ — Eji) 



771 717 

9l 9l '?77' 



For the terms with the factor t^e , only one is nonzero, that is 



2! 



since 



(iV'^-i*;D) = A2^'(iV'^v*;D) =0, 
AJ*' (i2e-i^-*;N) = AJ"' (i2e-'^-i*;D) = AJ'' (i^e-^^.'*; D) = 0. 



(A37) 



(A38) 

(A39) 

(A40) 
(A41) 



We can see that these terms can be absorbed into (or merged with) the lower order approximations to obtain the 
improved forms of perturbed solutions. 

2. 1=5 case 

Now let we consider the case of the fifth order approximation [l = 5). From ea.ip^ it follows that the first 
decompositions of ij-product have 2'^ = 16 terms. They can be divided into 5 groups 



(A42) 



i=0 
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where i indicates the number of rj functions. Obviously 

A2'''{0;v)^Ar'{cccc), 

-^V' i^iV) ^ AJ'^' (cccn) + AJ^' (ccnc) + AJ'^' (cncc) + A'p' (nccc) , 



AT (2;r/) 



AP (cnnc) 



AJ"' {ccnn) + AJ'^ (cncn) -r ^5 
+Aj'^'{nccn) + AJ'^' (ncnc) + A^'^'(nncc), 



/ / / / / 

^5''' (3; 77) = AJ'^ (cnnn) + AJ'^ (ncnn) + AJ'^ (nncn) + AJ'^ (nnnc), 

AJ'^ (4; 77) = AJ'^ {nnnn). 



(A43) 
(A44) 

(A45) 
(A46) 
(A47) 



Here, we have used the notations stated in Sec. I VII 

By calculation, we obtain the AP (0, rj) and every term of A^ (1, rj) have only nontrivial first contractions and/or 
anti-contractions. But, we can find that every term of AP (2, 77) can have one nontrivial second or third or fourth 
contraction or anti-contraction, that is 



477 


{ccnn) 


= AV{ 


AV 


{cncn) 


= Ar'{ 


AT 


{cnnc) 


= AT'i 


AT 


{nccn) 


= Ar'{ 


AT 


{ncnc) 


= Ar'{ 


AT 


{nncc) 


= Ar'{ 



177 



177 
^5 



{ncnc, ck) + A^ {ncnc, nk), 
c, ckk) + Ap {nncc, nkk). 

Similarly, every term of ^5'' (8,77) can have two higher order contractions and/or anti-contractions: 



(A48) 
(A49) 
(A50) 
(A51) 
(A52) 
(A53) 



Ap {cnnn) 
Ap {ncnn) 
Ap {nncn) 
Ap {nnnc) 



Ap {cnnn, kcc) -I- {cnnn, ken) 
+Ap {cnnn, knc) -\- Ap {cnnn, knn), 
Ap {ncnn, kkc, ck) -t- Ap {ncnn, kkn, ck) 
+Ap {ncnn, kkc, nk) + Ap {ncnn, kkn, nk), 
Ap {nncn, ckk, kc) + Ap {nncn, ckk, kn) 
+A2'^ {nncn, nkk, kc) + Ap {nncn, nkk, kn), 
Ap {nnnc, cck) + Ap {nnnc, cnk) 
+Ap {nnnc, nek) + AJ^ {nnnc, nnk). 



(A54) 
(A55) 
(A56) 
(A57) 



Moreover, their last terms, with two higher order anti-contractions, can have one nontrivial more higher contraction 
or anti-contraction: 

AJ^ {cnnn, knn) — Ap {cnnn, knn, kc) + Ap {cnnn, knn. kn), (ASS) 

Ap {ncnn, kkn, nk) ~ Ap {ncnn, kkn, nk,c) -\- Ap {ncnn, kkn, nk,n), (A59) 

Ap {nncn, nkk, kn) — AJ'^ {nncn, nkk, kn,c) + Ap {nncn, nkk, kn,n), (A60) 

Ap {nnnc, nnk) — Ap {nnnc, nnk, ck) + Ap {nnnc, nnk, nk). (-A-61) 

In the case of Ap {nnnn), there are three terms corresponding to the second decompositions that result in 

Ap {nnnn) = Ap {nnnn, ccc) -f Ap {nnnn, ccn) + Ap {nnnn, cnc) + A^'^ {nnnn, ncc) 

+A2'^ {nnnn, cnn) + Ap {nnnn, ncn) + Ap {nnnn, nnc) + Ap {nnnn, nnn). (A62) 
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In the above expression, from the fifth term to the seventh term have the third- or fourth- contraction and anti- 
contraction, the eighth term has two third contractions and anti-contractions: 



AJ'^ (nnnn, cnn) 
AJ'^ (nnnn, ncn) 
AJ"' (nnnn, nnc) 
{nnnn, nnn) 



AJ'^ {nnnn, cnn, kc) + A'P {nnnn, cnn, kn), 
Aip {nnnn, ncn, c) -f Al^ {nnnn, ncn, n), 
j4g^ {nnnn, nnc, ck) + AJ"' {nnnn, nnc, nk), 
{nnnn, nnn, cc) + AJ"' {nnnn, nnn, cn) 
+A2'^ {nnnn, nnn, nc) -f AJ'^ {nnnn, nnn, nn). 



In addition, AT^^ {nnnn, nnn, nn) consists of the fourth contraction and anti-contraction 

{nnnn, nnn, nn) — AJ'^ {nnnn, nnn, nn, c) -|- A'^'^ {nnnn, nnn, nn, n). 



(A63) 
(A64) 
(A65) 

(A66) 
(A67) 



According to the above analysis, we obtain that the contribution from the five order approximation is made of 52 
terms after finding out all of contractions and anti-contractions. 
Just like we have done in the / = 4 case, we decompose 



A, 



77 



Ar {e)+Ar {te)+Ar {t'e). 



where 



AT (e) 



AT (e-'^-'*) + AJ'' (e-'^-i*) + AJ^ (e-'^-^*) + AJ'' { 



Aj^'{te) = A^^'(te-'^-'*) + A^^'(te-'^-'i*) + A^^'(te-'^-'2t) + 
AJ^'{t^e) = ^^^'(t^e-'^-*) + A^T'(i2e-i^^-i*) + A^T'(t2e-i^^^2*) + A^V(^2g-i£;yt)^ 

While, every term in the above equations has its diagonal and off-diagonal parts about 7 and 7', that is 

^77'(g-iis,,t) ^ A^^'(e-'^^.*;D)+A^^'(e"'^-^*;N), 
A^^'(te-'^-.*) = Aj'''{te-'^-'^*;D)+Aj'''{te~'^-"';N), 
^^^'(i'e-'^-.*) = A^^'(t2e-'^-»*;D) + A^^'(t2e-'^-.*;N). 



(A68) 

(A69) 
(A70) 
(A71) 

(A72) 
(A73) 
(A74) 



where E-y. takes , E-y-^ , E^^ and E^i . 

If we do not concern the improved forms of perturbed solution higher than the fourth order one, we only need to write 
down the second and third terms in cq. (|A68p . We can calculate them and the results are put in the supplementary 
ofRef. [13. 

Based on these contraction- and anti contraction- expressions, we can, in terms of the rearrangement and summation, 
obtain 



^1 1^7 ~ -^^71 ) 



771 717 X , 
2i/l i/l "77 



71,72 



{E^ — Ey^ ) {Ej — Ey 



{E^ - E^J {Ey - E-f^) 



^771^7172^727^^^, ^ (-iGf 



(A75) 



where 



E 



71,72,73,74 

- E 

71,72,73 



„77l„7l72 7273 7374 747„ „ 
t/l Hi Hi ill Hi '/772'/773 

{Ey ~ EjJ {E^ - E^^) {E^ - E^^) {Ey - E^J 



„77l „772 „7l7 „7273 „737 
i/l i/l i/l i/l i/l 



„77l „772 „7l7 „72 73 ^737 

i/l i/l yi yi yi 



{Ey - Eyj'^ {Ey - E^^) {Ey - Ey^) {Ey - E^J {E^ ~ S^J^ {E^ - E^^) 



„77l „772 „7l7 „72 73 „73 7 

i/l yi iyi i/l yi 



{E^ - Ey^) {E^ - E^^) {E^ - £'^3) 



(A76) 
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71 (^7~-^7i) 



77' 



71,72 (-^'^ ~ -^Ti^ ^^■y-' ~ -^•^2) 



71,72 (-^7 ~ ^72) (-^71 ~ -^72) 



„77l „7l72 „727X 
ai i/l i/l "7 



(-^7 - -^7') 
-(-iGft)E 

71 L 



71 

g-i£,t^77l,^7l7^^77' 



^.-ii;-,tg77ig7l7'^77' 
(iJ^. — -^71) (-^7 ~ -^7') 



^-ii?-,t^^77i ^^717^^77 



(^/^ — -^71) (^7 — -^7') (-^7 — -^71) (-^7 — -^7') 
e-ii^.*^77l^7l72^727'^^^^^^^, 



7i>72 



(-£^7 ~ -^71 ) (-^7 ~ -^72 ) (-^7 ~ -^7' ) 



(A77) 

(A78) 
(A79) 



(A80) 



A5(te-'^-*,N) = -(-it)E 

71 

-(-io E 



4^^-'^^i*ff7^^ff7^^'r?-,y 
{E^ - E^J {E^^ - Ey) 

{E^ - Ej^) {Ej^ - E^^) {E^^ - Ej>) 



A5{te-'^y-*,N) = i-it) E 

71,72 



[E-i - -E72) (-^71 ~ -^72) (-^72 - ^i') 



71 



{Ej — -^7') (-^71 — -^7') 



+(-iG?^i)E 



71 



^^^0 E 



/t 7'7i„7i7' 77' o-i-^-v'*n'l"7i „7i7'„77' 

^ ^ ffl ffl ffl |_ ^ ^ ffl ffl ffl 

(£/7 — -^7') (-^71 — -^7') (-^7 — -^7') (-^71 — -^7') 

,-i^,,tg77ig7i72g727'^^^^,^^^, 

- -E7') (£^71 - -B7') (-B72 - -E7') 



For the parts with f^e, we have 



A5(t2e-'-^^*,D) 
A5(t^e-'^^i*,D) = ^5(t2e-'^^2*,D) =^5(t2e-'^v*,D) =0. 



2! 



■gV 



A5(t2e-'^^2*,N) = A(i2c-'^v*,N) =0 

2! V T' ; p _ p , 



A5(i'e-'-^v*,N) 



,2 e-'^7'* 



57^ 



(A81) 
(A82) 



(A83) 

(A84) 
(A85) 

(A86) 
(A87) 
(A88) 



It is clear that the above diagonal and off-diagonal part about {te) and AJ'^ (te) indeed has the expected 
forms and can be absorbed reasonably into the lower order approximations in order to obtain the improved forms of 
perturbed solutions. 
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3. 1 — 6 case 



Now let we consider the case of the sixth order approximation {I — 6). From eq. flS]) it follows that the first 
decompositions of g-product have 2'^ = 32 terms. Like the I = 5 case, they can be divided into 6 groups 

^Y^^r'i^iv), (A89) 

i=0 

where i indicates the number of rj functions. Obviously 

A'^/ (0 ; r;) = A^'^' (ccccc) , ( A90) 



AJ"^ (1;??) = AJ'^ (ccccn) + AJ'' {cccnc) + AJ'^ {ccncc) 

+AJ'^' (cnccc) + Al''' {ncccc) , (A91) 



AJ'^ (2; r/) = AJ'' {cccn) + AJ"^ {ccncn) + AJ'' {cnccn) + AJ'' {ncccn) 

+Aq'^ (ccnnc) + A'^'^ [cncnc) + A'^'^ (nccnc) + AJ^ (cnncc) 

+Al'^' [ncncc) + AJ^' {nnccc) , (A92) 



AJ^ (3; 77) — Aq'^ (ccnnn) + (cncnn) + AJ^ (cnncn) + Aq'^ (cnnnc) 
+ylg^ (nccnn) + AJ^ (ncncn) + A^'^ (ncnnc) + AJ^ (nnccn) 

+A2'^' (nncnc) + AJ'^' innncc) , (A93) 



Ap (4; if) — A^'^ (cnnnn) + AJ'^ (ncnnn) + A'^ (nncnn) 

+Ag^ {nnncn) + Aq'^ (nnnnc) (A94) 



AJ"'' (5; 77) = AJ"'' (nnnnn) . (A95) 
Furthermore considering the high order contraction or anti-contraction, we have 

Ap (cccnn) = Ap (cccnn, kkkc) + A^'^ {cccnn, kkkn), (A96) 

Ap (ccncn) — Ap (ccncn, kkc) + A^ (ccncn, kkn), (A97) 

Ap (ccnnc) = Ap (ccnnc, kkck) + Ap (ccnnc, kknk), (A98) 

/ / / 

Ap (cnccn) — Ap (cnccn, kc) + Ap (cnccn, kn), (A99) 

Ap (cncnc) — Ap (cncnc, kck) + Ap (cncnc, knk), (AlOO) 

Ap (cnncc) — Ap (cnncc, kckk) + AJ'^ (cnncc, knkk), (AlOl) 

Ap (ncccn) — Ap (ncccn, c) + Ap (ncccn, n), (A102) 

Ap (nccnc) = Ap (nccnc,ck) + Ap (nccnc,nk), (A103) 

Ap (ncncc) = Ap (ncncc, ckk) -'r Ap (ncncc, nkk), (A104) 

Ap (nnccc) = Ap (nnccc, ckkk) + Ag^ (nnccc, nkkk) (A105) 



Ap (ccnnn) = Ap (ccnnn, kkcc) + A^'' (ccnnn, kkcn) + AJ'^ (ccnnn, kknc) 

-'rAp (ccnnn, kknn, kkc) + Ap (ccnnn, kknn, kkn), (A106) 
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Aq'^ {cncnn 

I 

JC^ {cnncn 
Aq'^ {cnnnc 

A'P (nccnn 
Aq'^ {ncncn 
Aq'^ {ncnnc 
Aq'^ {nnccn 
Aq'^ {nncnc 
AJ"^ {nnncc 



Aq'^ {cncnn, kkkc, kck) + AJ^ {cncnn, kkkc, knk) + Al'^ {cncnn, kkkn, kck) 

{cncnn, kkkn, knk, kc) + A^^ {cncnn, kkkn, knk, kn), 
A'ip {cnncn, kckk, kkc) + AJ^ {cnncn, kckk, kkn) + AJ^ {cnncn, knkk, kkc) 
+AJ'^ {cnncn, knkk, kkn, kc) + A^"' {cnncn, knkk, kkn, kn), 
A'P {cnnnc, keck) + A'^ {cnnnc, kcnk) + AJ'^ {cnnnc, knck) 
+Aq'^ {cnnnc, knnk, kck) + AJ^ {cnnnc, knnk, knk), 

AJ^ {nccnn, kkkc, ck) + AJ'^ {nccnn, kkkc, nk) + AJ^ {nccnn, kkkn, ck) 

{nccnn, kkkn, nk, c) + AJ"' {nccnn, kkkn, nk, n), 
Aq ' {ncncn, ckc) + A^ ' {ncncn, ckn) + {ncncn, nkc) 
+AJ'^ {ncncn, nkn, c) + AJ'^ [ncncn, nkn, n), 

Aq'' {ncnnc, kkck, ckk) + AJ'' {ncnnc, kkck, nkk) + AJ'' {ncnnc, kknk, ckk) 
+Aq^ {ncnnc, kknk, nkk, ck) + Ap {ncnnc, kknk, nkk, nk), 
A'P {nnccn, cc) + A'P {nnccn, cn) + Af {nnccn, nc) 
+Aq"' {nnccn, nn, c) + A'P {nnccn, nn, n), 

Aq^ {nncnc, ckkk, kck) + A^ {nncnc, ckkk, knk) + Aq'^ {nncnc, nkkk, kck) 
+Aq'^ {nncnc, nkkk, knk, ck) + AJ'^ {nncnc, nkkk, knk, nk), 
A'p {nnncc, cckk) + AJ'^ {nnncc, cnkk) + AJ'^ {nnncc, nckk) 
+Aq^ {nnncc, nnkk, c) + Ap {nnncc, nnkk, n) 



(A107) 
(A108) 
(A109) 

(Alio) 
(Alll) 
(A112) 
(A113) 
(A114) 
(A115) 



Aq'^ {cnnnn) = A^ {cnnnn, kccc) + A^ {cnnnn, keen) + Ap {cnnnn, kcnc) 
lit 
+Ap {cnnnn, knee) + A^'^ {cnnnn, kenn) + Ap {cnnnn, kncn) 

I I 

-\-Ap {cnnnn, knnc) + Ap {cnnnn, knnn). 



Ap {cnnnn, kenn) = Ap {cnnnn, kenn, kkc) + Ap {cnnnn, kenn, kkn). 



Ap {cnnnn, kncn, kc) + Ap {cnnnn, knen, kn), 
Ap {cnnnn, knnc) = Ap {cnnnn, knnc, kck) + Ap {cnnnn, knnc, knk) 



Ap {cnnnn, kncn) 
177' 



(A116) 

(A117) 
(A118) 
(A119) 



Ap {cnnnn, knnn) = Ap {cnnnn, knnn, kcc) + Ap {cnnnn, knnn, ken) 

I I 

+Ap {cnnnn, knnn, knc) + Ap {cnnnn, knnn, knn, kc) 
+Ap {cnnnn, knnn, knn, kn). 

Ap {ncnnn) = Ap {ncnnn, kkcc, ckk) + Ap {ncnnn, kkcc, nkk) 
+Ap {ncnnn, kkcn, ckk) + Ap {ncnnn, kknc, ckk) 
+Ap {ncnnn, kkcn, nkk) + Ap {ncnnn, kknc, nkk) 
+Ap {ncnnn, kknn, ckk) + Ap {ncnnn, kknn, nkk). 



Ap {ncnnn, kkcn, nkk) 

Ap {ncnnn, kknc, nkk) 
I 

A^'^ {ncnnn, kknn, ckk) 



Ap {ncnnn, kkcn, nkk, c) + Ap {ncnnn, kkcn, nkk, n), 
Aq ' {ncnnn, kknc, nkk, ck) + Ap {ncnnn, kknc, nkk, nk). 



(A120) 



A^'^ {ncnnn, kknn, ckc) 



Ap {ncnnn, kknn, ckn). 



(A121) 

(A122) 
(A123) 
(A124) 



{ncnnn, kknn, nkk) = A^'^ {ncnnn, kknn, nkc, ck) + AJ'^ {ncnnn, kknn, nkc, nk) 

+AJ'^ {ncnnn, kknn, nkn, ck) + A'q' {ncnnn, kknn, nkn, nk, c) 
+Ag''' {ncnnn, kknn, nkn, nk, n) 



A'l'^ {nncnn) 



AJ^ {nncnn, ckkc, kck) + A'P {nncnn, ckkc, knk) 
+Ag^ {nncnn, ckkn, kck) + {nncnn, nkkc, kck) 
+A1'^ {nncnn, ckkn, knk) + Ag'^ {nncnn, nkkc, knk) 
^-Aq'^ {nncnn, nkkn, kck) + AJ^ {nncnn, nkkn, knk). 



Aq'' {nncnn, ckkn, knk) 
Aq'^ {nncnn, nkkc, knk) 
A'P {nncnn, nkkn, kck) 



Aj'^ {nncnn, ckkn, knk, kc) + Af {nncnn, ckkn, knk, kn), 
AJ'^ {nncnn, nkkc, knk, ck) + Aq'^ {nncnn, nkkc, knk, nk), 
A'P {nncnn, nkkn, kck, c) + AJ'^ {nncnn, nkkn, kck, n). 



AJ'^ {nncnn, nkkn, knk) = AJ'' {nncnn, nkkn, knk, cc) + AJ'' {nncnn, nkkn. knk, cn) 

+A2^ {nncnn, nkkn, knk, nc) + ^g'*' {nncnn, nkkn, knk, nn, c) 
+Ag''' {nncnn, nkkn, knk, nn, n) 



A'l'' {nnncn) 



Al'' {nnncn, cckk, kkc) + A^'' {nnncn, cckk, kkn) 
+AJ^ {nnncn, cnkk, kkc) + AJ^ {nnncn, nckk, kkc) 
+A'P {nnncn, cnkk, kkn) + Af^ {nnncn, nckk, kkn) 
+A2^ {nnncn, nnkk, kkc) + AJ^ {nnncn, nnkk, kkn). 



AJ'^ {nnncn, cnkk, kkn) = AJ'^ {nnncn, cnkk, kkn, kc) + A'P {nnncn, cnkk, kkn, kn), 

Aq'^ {nnncn, nckk, kkn, c) + AJ'^ {nnncn, nckk, kkn, n), 
A'P {nnncn, nnkk, ckc) + AJ"* {nnncn, nnkk, nkc). 



Ag'^ {nnncn, nckk, kkn) 
A'P {nnncn, nnkk, kkc) 



Aq {nnncn, nnkk, kkn) 



AJ^ {nnncn, nnkk, ckn, kc) + AJ'^ {nnncn, nnkk, ckn, kn) 
+AJ^ {nnncn, nnkk, nkn, kc) + AJ^ {nnncn, nnkk, nkn, kn, c) 
+A'q* {nnncn, nnkk, nkn, kn, n). 



Ag'*' {nnnnc) = Ag'*' {nnnnc, ccck) + A'f {nnnnc, ccnk) + Aj'' {nnnnc, cnck) 

+A'P {nnnnc, neck) + AJ'^ {nnnnc, cnnk) + AJ'^ {nnnnc, ncnk) 
+A'P {nnnnc, nnck) + AJ'^ {nnnnc, nnnk), 

Aq'^ {nnnnc, cnnk) = Aq'^ {nnnnc, cnnk, kck) + AJ'^ {nnnnc, cnnk, knk), 
A'P {nnnnc, ncnk) = AJ'^ {nnnnc, ncnk, ck) + AJ'^ {nnnnc, ncnk, nk), 
A'P {nnnnc, nnck) = AJ'^ {nnnnc, nnck, c) + AJ'^ {nnnnc, nnck, n), 



{nnnnc, nnnk) 



A'P {nnnnc, nnnk, cck) + {nnnnc, nnnk, cnk) 
+A'P {nnnnc, nnnk, nek) + A'P {nnnnc, nnnk, nnk, ck) 
+Aq^ {nnnnc, nnnk, nnk, nk). 
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A'g' (nnnnn) = 



Ap {nnnnn, cccc) + Ap {nnnnn, cccn) + {nnnnn, ccnc) 

/ / f 

+A^'^ [nnnnn, cncc) + Ap 



, ^ , ^ , ^ , ^ , 

[nnnnn, cncc) + Ap [nnnnn, nccc) + A^ [nnnnn, 
-\-Ap [nnnnn, cncn) + A^ [nnnnn, cnnc) + [nnnnn 
+Ap [nnnnn, ncnc) + A^ [nnnnn, nncc) + A^f (nnnnn 



ccnn) 
nccn) 



cnnc) + Ap [nnnnn, 

cnnn 

+Ap [nnnnn, ncnn) + A^ [nnnnn 
^-Ap [nnnnn, nnnn). 



j -r Ap [nnnnn, 
[nnnnn, ncnn) + Ap [nnnnn, nncn) + Ap [nnnnn, nnnc) 



Ap [nnnnn, ccnn 
Ap [nnnnn, cncn 
Ap [nnnnn, cnnc 
Ap [nnnnn, nccn 
Ap [nnnnn, ncnc 
Ap [nnnnn, nncc 



Ap [nnnnn, ccnn, kkc) + Ap [nnnnn, ccnn, kkn), 
Ap [nnnnn, cncn, kc) + Ap [nnnnn, cncn, kn), 
Ap [nnnnn, cnnc, kck) + Ap [nnnnn, cnnc, knk), 
Ap [nnnnn, nccn, c) + Ap [nnnnn, nccn, n), 
Ap [nnnnn, ncnc, ck) + Ap [nnnnn, ncnc, nk), 
AJ'^ [nnnnn, nncc, ckk) + Ap [nnnnn, nncc, nkk). 



(A141) 

(A142) 
(A143) 
(A144) 
(A145) 
(A146) 
(A147) 



Ap [nnnnn, cnnn) = Ap [nnnnn, cnnn, kcc) + Ap [nnnnn, cnnn, ken) 

+Aq'^ [nnnnn, cnnn, knc) + Ap [nnnnn, cnnn, knn, kc) 
+A2'^ [nnnnn, cnnn, knn, kn), 
Ap [nnnnn, ncnn) = Ap [nnnnn, ncnn, kkc, ck) + Ap [nnnnn, ncnn, kkc, nk) 

■^Ap [nnnnn, ncnn, kkn, ck) + Ap [nnnnn, ncnn, kkn, nk, c) 
+Ap [nnnnn, ncnn, kkn, nk, n), 
Ap [nnnnn, nncn) = Ap [nnnnn, nncn, ckk, kc) + Ap [nnnnn, nncn, ckk, kn) 

+Ap [nnnnn, nncn, nkk, kc) + Ap [nnnnn, nncn, nkk, kn, c) 
+j4g^ [nnnnn, nncn, nkk, kn,n), 
Ap [nnnnn, nnnc, cck) + Ap [nnnnn, nnnc, cnk) 
+Ap [nnnnn, nnnc, nek) + Ap [nnnnn, nnnc, nnk, ck) 
+Ap [nnnnn, nnnc, nnk, nk), 

Ap [nnnnn, nnnn, ccc) + Ap [nnnnn, nnnn, ccn) 
+Ap [nnnnn, nnnn, cnc) + Ap [nnnnn, nnnn, ncc) 
+Ap [nnnnn, nnnn, cnn) + Ap [nnnnn, nnnn, ncn) 
+Ap [nnnnn, nnnn, nnc) + Ap [nnnnn, nnnn, nnn). 



Ap [nnnnn, nnnc) 



Ap [nnnnn, nnnn) 



(A148) 



(A149) 



(A150) 



(A151) 



(A152) 



Ap [nnnnn, nnnn, cnn) = Ap [nnnnn, nnnn, cnn, kc) + Ap [nnnnn, nnnn, cnn, kn), 
Ap [nnnnn, nnnn, ncn) = Ap [nnnnn, nnnn, ncn, c) + Ap [nnnnn, nnnn, ncn, n), 
Ap [nnnnn, nnnn, nnc) = Ap [nnnnn, nnnn, nnc, ck) + Ap [nnnnn, nnnn, nnc, nk). 



(A153) 
(A154) 
(A155) 



Ap [nnnnn, nnnn, nnn) = Ap [nnnnn, nnnn, nnn, cc) + Ap [nnnnn, nnnn, nnn, cn) 

I I 

+Ap [nnnnn, nnnn, nnn, nc) + Ap [nnnnn, nnnn, nnn, nn, c) 
+Ap [nnnnn, nnnn, nnn, nn,n). 



(A156) 



Thus, we obtain that the contribution from the six order approximation is made of 203 terms after finding out all 
of contractions and anti-contractions. 
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Just like we have done in the ^ = 4 or 5 cases, we decompose 

AT' = Al''{e) + Al<'(te) + A}''' {fe) + Al<' {t'e) (A157) 
= Al''{e,te)+Al'\ee,r'e), (A158) 

To our purpose, we only calculate the second term A'l'^ (t^e, t^e) in ea. (jA158[) . Without loss of generality, we decompose 
it into 

While, every term in the above equations has its diagonal and off-diagonal parts about 7 and 7', that is 

AJ'' {t'e-''^-'^') = A2^'(i2e-'^-.*;D)+A;r^'(i2e-'^-.*;N), (A160) 

Al'\t^e-''^-^') = Al<' {t^e-''^-^^';V)) + Al<' {t^e-''^-^'-^). (A161) 

where E~f^ takes Ey,E-y-^ and E^i. 

Based on our calculations, we find that there are nonvanishing 91 terms and vanishing 112 terms with t'^e, t^e factor 
parts in all of 203 contraction- and anti contraction- expressions (see in the supplementary Ref. [13]). Therefore we 
can, in terms of rearrangement and summation, obtain the following concise forms: 



2 



(-it)- (g?^ 

--^y - — ^gV'gT^^^i'- (A162) 

2! ^ {E,-E,,f 

^ =.^2 fc^'^Ve-i^-i* 
A6(t^e-'^-*,D) = 4f^E \^ ^ ^ ,2 g7^^grV- (A163) 

A(i^e-'^-''*,D) = 0. (A164) 
AJf'-e-'^-'*^ n\ - i:il^2(7(2)(7(3)_£j^Z^„77' 

^6(i e ,N) - ^2G^, G^, ^E^_E^, fi 

Their forms are indeed the same as expected and can be absorbed reasonably into the lower order approximations in 
order to obtain the improved forms of perturbed solutions. 
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Supplement 



11 — 5 case 

In the following, we respectively calculate the 52 component expressions of AJ'^ , and put the second and third 
terms in ea. (|A68p together as AJ'^ (te, i^e). 



{cccc; te, t e) = 



-(-it) 



3e 



3e 



{Ej - Ey) 



2 {E,-Eyf 



iE^jt 



2 {E,-Eyf 



11 



11 



(168) 



A'P {cccn;te,t^e) = ^ 



-(-it) 



2e- 



AE-yt 



{Ej — E^-^ ) {E^ — E^r ) 



(-it) 



-iE~,t 



-(-it) 



jj -^71 



{E^ — ^71 ) (-^7 — -^7') 



(-it)2 



71/ 
iE^t 



e ^ 



{E-y - E^,y {E^-, - Ey) 2! {E-y - E^,y [E^ - Ey) 



X \gr'\ gV »77i7'- 



(169) 



AJ* {ccnc; 



(-it) 



fyEy — -E-yi ) {E^ — E^i'y' (-^7 — -^7') (-^71 ~ -^7') 



-(-it) 



2e 



-iE^it 



^-iEyt 



{E-y — Ej') (-^71 -^1 



2! 



/ 2 ,2 , 

X gf^ gV gV v^ji- 



(170) 



AJ'' {cncc; te, t^e) — ^ 



71 L 



-(-it) 



2e-' 



E-y t 



i_E_,t 



{Ey — £/-yi ) {Ej — E^i ) 



(-it) 



Jy -^71 

e-'^v* 



(-it)2 



a-iE^t 



(Ey — E^i) {Ej-^ — ^7') ^' (-^7 ~ -^71 ) (-^7 — ^7') 



X \g 



771 I 



gi 



77 



77 

5l ^771 7' 



(171) 



37 



AJ"' {nccc;te,t^e) = 



71 L 



-i-it) 



2e- 



3 - {-it) 3 

{Ej — -B-yi ) (-£'71 — -£'7') (-£'7 — E-y') ~ -£'7') 



-i-it) 



2! 



(^//y -^7') (-^71 -^7') 



^717 



51 ^771 • 



(172) 



{ccnn, kkc; te, t^e) 



E 



7i>72 L 



-(-ii) 



-(-it) ^-"2 2 ~ (~^^) ~~3 

{Ef — -E71) {E^ — E^^) {Ey — E^J {E^ - E^^) 



+ 



(^7~-^7l) (^71^-^72) ^' i^l - E^l) (-£^7 --^72) 



x|5n'<?7''^5?'''5r^'577'- 



(173) 



{ccnn,kkn;te,t'^e) = 



71,72 



(-it) 



{E^ - E^J {Ej - E^^) {Ej - E^>) 



+ (-«) 



(-£^7 ~ -^71)^ (-^71 ~ -^72) (■£'71 ~ -^7') 



X bi^'l gV^gT^^gT^ %i7'??772??77'- 



(174) 



{cncn,kc:tc,t'^c) 



E 



( 'i) 2c"'^^* ^ .^^ e~'^i'* 

(ii/7 — £"71) (-£'7 — Ej/)"^ {Ej — -£'71) (-£'7 — -£^7') 



+ (-ii) 



x\gT 



+ 



(-ii)^ 



{E^ — -^7') (-^71 — -^7') ^' (-^7 — -^71) (-^7 — -^7' 
2 



77 

gr 



77 

5l ^717' 



(175) 



AJ"^ {cncn,kn;te,t^e) = 



-(-it) 



Q—iE^t 



7i>72 _ 

-(-it) 
-(-it)- 



{E^ — E^^) {Ej — E^^) {Ej - Ey) 
{E^ - E^J [E.f - E^i^f {E^ - Ey) 



+ 



{E-y - J^7i ) {Ej - E^^ ) {Ej - Ey ) 



2! {E^ — Ej^) {Ej — E^^) {Ej - Ey) 



I 5i ^7172 ^71 7' ^72 7' 



(176) 



Ag {cnnc, kck;te,t e) 



= E 



Q—iE^t ^—iE^-^t 

-it) 2 2 + 2 2 

(£/7 — -£^71) (-£'7 — -£'7') (-£^7 — -£'71) (-£'71 — -£'7') 

2 



+(-ii) 



{E.y — E^l) (-E71 -£'7') 



\gj 



111 I 



^717 



gV 



(177) 
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{cnnc, knk;te,t^e) = 



7l>72 L 



i-it) 



{E^ - E^^ ) (Ej - Ej^ ) {Ej - Ey ) 



+{-it) 2- 

{E^ - Ej>) {E^^ - Ery,) - E^i) 



771 1 2 



^7.7 



Si ^7l7''?7l72'?772- 



(178) 



Aj"' {nccn, c;te,t^e) 



E 



(-it) 



(-E'7 ~ ^71)^ (-^7 ~ -^72)^ ^ ' ^ (-^7 ~ -^71)^ (-^71 ~ -^72)^ 



(-^7 ^ ^72)^ (-^71 ~ -^72)^ 



1^7172 |2 ^771^7172^727^^^,. 



(179) 



{nccn,n;te,t^e) = 



-(-it) 



(-E'7 ~ -^71) (-^71 ~ -^72) (-S71 ~ -^7') 



-(-it) 

(-^7 ~ -^72) i^-fi - -^72)^ (-£^72 - -E'7') 

^ l^i^'^^l 9i^^ 9i^^^ 9i^^ ^71 7' ^772 ^77'- 



(180) 



Ag''' {ncnc,ck;te,t^e) 



= E 



71 



(-it) 



^ .^^ e-'^7'* 

(£/-y -^71) (-^7 -^7')^ (-^7 -^7') (-^71 -^7') 



-(-it) 



(i/zy — E^f^^ (-^71 — -^7') ^* (-^7 — -^7') (-^71 — -^7') 



7i7 

51 



77 

5i 



5i ^771 • 



(181) 



(ncnc,nk;te,i^e) — 



7i>72 



-it) 



-(-it) 



-(-it) 



(-£^7 ~ -^7') (-^71 ~ -^7') (-^72 ~ -^7') 

(-£^7 - -E'7') (-E71 - -E^')^ (-B-y2 ~ -^7') 



(i;^ - Ey) {E^, - Ey) {E.,, - E.,,) 



(-it)2 



2! {Ej - Eji) {E^^ - Ey) (£72 ~ -^7') 
/ 2 ,2 , 
X g';^ g-^ 'ri"cyi'njj2Vji-y2- 



(182) 



A5''' {nncc,ckk;te,t^e) 



= E 

7ii72 



-(-it) 



2e 



-iE-,t 



{Ej - EjJ {Ej - E^^) 



-lE-yt 



{Ej - E^^f {E^ - Ej^f 



(-it)^ 



+(-it) 3 ■ 91 2 

{E^ - EyJ {Ej^ - Ej^) ^- {E^ - E^J {Ej - Ej^) 



(183) 
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{nncc,nkk;te,t^e) = 



i-it) 



71)72 L 



(-^7 ~ -^72) (-^71 ~ -^72) (-S72 - ^1') 



(-E'7 ~ -^7') (-^71 ~ -^7') (-^72 ~ -^7') 



^ Si^''' g'T'^ gl^'^^ gl^'^ ^77i7''?772 ^77' 



(184) 



^5''' {cnnn, kcc; te, t^e) 



E 

7i>72 



-(-it) 



-(-«) 



{Ej - Ej^) {E^ - Ej^) 



(-£^7 ~ ^lif (-^7 ~ -^72)^ 



2 



2! 



(-^7~-^7i) i^l - ^12) 



x\grrgrgr'grs,^'- 



(185) 



A5''' {cnnn, ken; te,t^e) = ^ 



(-it) 



{Ej - EjJ {E^ ~ E^^) {Ej - Ejr) 

(^7 - Eflf (^71 ~ -^72) (-^71 ~ -^7') 



^|^77.|^^772^727l^7.7^^^^,^^^,. 



(186) 



[cnnn, knc; te, t^e) = ^ 

7l>72,73 
-(-it) 

-(-it)- 



-(-it) 



(-E'7 - -2^71 ) (-^7 " -^72 ) {^1 - -E'73 ) 

^—'lEyt 



{Ej - Ej^ ) {E^ - Ey^ ) {E^ - Ej^ ) 



+ 



{E-Y - -Et,i) (-E^ - E^^) {E^ - E^^) 
{-itf e-'^^* 



2! (-E'7 — -£'71) (-£"7 ^ -^72) (-^7 ^ -^7 



^ I fll '''^ S'l^^'771 72 '?7l 73*^77' ■ 



(187) 



{cnnn, knn, kc; te, t^e) = 



71.72 



(-it) 



{Ej - E^^ ) (£^ - E^^ ) (E^ - Ey 
{Ej - Ej')^ {Ej^ - Ey) (£^2 - £y ) 



/ 2 , 

X gV ??772%i7' 



(188) 



aI ' {cnnn, knn, kn; te, t e) 



7i>72,73 



e ' Iffl'*' I gT^ gj gj ^773 ^77'^7l72^7l73^7l7'^^727' 

{E^ - E^J {E^ - E^^) {E^ - E^^) {E^ - E^) 



(189) 



40 




41 



{nncn, nkk, kc; te, t^e) = 



7ii72 



-(-it) 



{E^ - E^J {E^^ - E^J {E^^ - Ey) 



+(-ii) 



{E^ - Ej') {E^^ - Ey)^ {E^2 - Ey) 



(197) 



Ap {nncn, nkk, kn, c; te, t^e) = 

->72' 

-(-it) 



7i>7273 L 



(-it) 



Q—iE^t 



{E^ - EjJ{E^- E^^ ) {E^- E^^) 



{E^ - E^^f {E^^ - E^^) {E^^ - E^^) 



(198) 



AJ'^ {nncn, nkk, kn, n; te, t^e) 

^ I 'y273 |2 771 7172 727' 

_ ^ -^-^ S ^2 I ^772%73'^77''?7l73^7l7'^737' 

^.■r:^^, (-^7 " -^72) (-E'71 ~ -^72) (-^72 - -^73) (-^72 - -^7') 



(199) 



{nnnc, cck; te, t e) 



E 



71)72 L 



-(-it) 



(E"^ - -E-yi)^ (E^ - E^a)^ (E^ - E^-i)^ (-Bt. - 



-(-it) 



^—'lEy^ t 



+ 



{E^ - EjJ {E^^ - E^^) 2! _ E^^) {E^ - E^^) 



X 1^771 |2 ^771 ^7172^727^^^,_ 



(200) 



AJ^ {nnnc,cnk;te,t^e) = 

71)72,73 

-(-it) 
-(-it)- 



-(-it) 



[E-y - E^J {E^ - Ej^) {E^ - Ej^) 

Q—iEyt 



{Ey - E^J [E^ - E^i^f {E^ - E^^) 



+ 



{E-Y - -E^i ) {Ej - -E^2 ) {E-Y - -E^a ) 

(-it)2 



Q—iEyt 



2! {E^ - E^^) {Ej - E^^) {Ej - Ej^) 

X \9i'^^\ 9i^'^^ 91^'^ 'niii3'ni2i3^ii' ■ 



(201) 



AJ^ {nnnc, nek; te, t^e) = 



7i>72 L 



-(-it) 



+(-ii) 



{Ej - E^J {Ej^ - Ey^) {E^^ - Ey) 

(Ef - -By) {E^^ — Eyf' {Ej^ - E^i) 
2 / 

„77i„7i72 727 „ „ , 
91 9l 9l '7772 '/77- 



(202) 



42 



{nnnc,nnk, ck;te,t^e) = 



7i>72 



i-it) 



{E^ - E^^) {Ey - E^^) {Ej - Eji) 



2- 

{Ef - -E7') (£'71 - E^>) {Ej2 - Ej') 



(203) 



AJ'^ {nnnc, nnk, nk; te, t^e) 



^737 



^77 gJ^g'J'^ ^772'7773^77'^7i73^7i7'^7273 



(-2^7 ~ -^7') (-^71 ~ -^7') (-^72 ~ E-r') (-^73 ~ -^7 



(204) 



A2'' {nnnn,ccc;te,t^e) 



{-it) — 2 + (~^^) 2 2 

(i/^ — -^71) (^7 — -^7') (-^7 — -^71) (-^71 — -^7') 



(i/^ — -^7') (-^71 — -^7') 



/\ 2 



1^771 |2 (^^771^717- j ^77' 



(205) 



Ag'*^ {nnnn,ccn;te,t'^e) = ^ 



(-ii) 



71,72 L 



{Ej-EjJ {E^ - E^^) {Ej - Ey) 



Q—iE^^ t 

{E^ - E^^f {E^^ - £^7 J (-5^71 ~ -^7') 



X(57^^5r^^ffr^)57^^5?"^77'^727'- 



(206) 



Ag''' {nnnn,cnc;te,t^e) = ^ 



i-it) 



7l>72 L 

( 



(£7 - £71) (£7 - £72) (£7 - £7') 
(£7 - £7') (£71 - £7') (£72 - E7/) 



XlflZ-'^ff?^"') (57^^57="') 57"'^7i72- 



(207) 



vlg''' {nnnn,ncc;te,t^e) = 



-(-it) 



(£7 - £7 J (£7^ - £"72) (£7^ - £7/) 



(£7 - £7') (£71 - £7') (£72 - £7') 

) 



(208) 



AJ'^ {nnnn,cnn,kc;te,t'^e) = 



(-it) 



(£7 - £71 ) (£7 - £72 ) (£7 - £7' ) 
(£7 - £7')^ (£71 - £7') (£72 - £7') 



x{9t'''9r) {9r9r'')9r'v,^,.- 



(209) 



A'^'^ (nnnn, cnn, kn; te, t^e) 



-iB-t „77l „7l72 _727 „773 „737 



• 5l 9l 9l 9l 9l '?77'^7i73^7273^7l7''7727' 



{E^ - E^^) {E^ - E^^) {E^ - E^^) {Ej - Ey) 



Ap {nnnn,ncn,c;te,t^e) = 

7l>7273 



i-it) 



Q—iEyt 



{E^ - Ej^f {E^ - Ej^) {E^ - Ej^) 



Q—iEy^ t 



AlI (nnnn, ncn, n; te, t e) 



71,72,73 



C ^1 gj"^ gj gj '^^gj-^ gj ^ ??77'^7772 ^773 ^727'^737' 

{Ey - EjJ {E^^ - £72) (£71 - E^^) {Ej^ - Ey) 



{nnnn, nnc, ck; te, t^e) = 



71,72 



i-it) 



+(-ii) 



{Ej - E^^ ) {E^ - E^^ ) {E^ - Ey ) 

{Ej - Ey) {Ej^ - Ey) {Ej^ - Ey) 



y-{9V'9r') (^r^^r) 57"'^7i72- 



Ap {nnnn, nnc, nk; te, t'^e) 



^ ^^^ ^ '' 9i ^ 9i ^ 9i^^ 9i^ 9i ^77^^772^773^7172^7173 



Ap {nnnn, nnn, cc; te, t^e) = 



71,72 



{-it) 



+ (-ii) 



{Ej - Ej^) {Ej — Ej^) {Ej — Eji) 
{Ej - Eji'f {Ej^ - Eji) {Ej^ - Eji) 



^gv 9r'9r'9r9r 



Ap {nnnn, nnn, cn; te, t^e) 



E (-i^) 



71,72,73 



e ^^yt gj'^"'^ gj^'^^ gj^"' gP ??772^7i73^7i7'^727'^737' 

{Ej - Ej^) {Ej - Ej^) {Ej " Ej^) {Ej - Ejl) 



Ap {nnnn, nnn, nc: te, t e) 

E -i-E /t 7'71„7172 7273 737' 77' nnnn, 
(^ — iiy ^1 ^1 ^1 ^1 -^1 '/771 '/772 '/773 '/7l73'/727 



{Ej - Ej>) {Ej^ - Ej>) {Ej2 - Ej,) {Ej^ - Ej,) 



Ap {nnnn, nnn, nn, c; te, t^e) 



E (-i*) 



71,72,73,74 



6 ffi 9i 9i 9i 9i ^772^773^7173^7174^7274^77' 

{Ej - Ej^) {Ej - Ej^) {Ej " Ej^) {Ej - EjJ 



Ap {nnnn, nnn, nn, n; te, t^e) = 



44 



D. 1 = 6 case 

Based on our calculations, we find that there are nonvanishing 91 terms and vanishing 112 terms with t^e, t^e factor 
parts in all of 203 contraction- and anti contraction- expressions. In the following, we respectively calculate them 
term by term, and we only write down the non-zero expressions for saving space. 



Ae{ccccc;t^e,t^e) = 



71 L 



e 



-iE~,t 



3! {E-y - 



(-It) 



2 Q—iEy^t 



1^ 2 {E^-E^,f 



(220) 



AQ{ccccn;t^e,t^e) = 



71 



Q—iEyt 



71/ V-^7 



e '-^71 



{E^ - EjJ (E^i - Ey) 



(221) 



Ae{cccnc; t^e,t^e 



71.72 L 



i-itf 



-lE^t 



3! {Ej - Ejj^) {E^ - E^i) 2 [E^ - E^^) {E^ - Ej^) 



-iE-,t 



2 {Ey - Ej^) {E-y - E^^) 



\9i 



""^''|57"1^717.^77'- 



(222) 



A%{ccncc\t^e,t'^e) = 



(-it)' 



Q—iE-^-^ t 



71 



2 {E^ - E^if {E^-L - Ey) 



gV'gr %7' 



(223) 



AQ{cnccc;t^e,t^e) = 



71.72 



-iE^t 



-iE^t 



3! {Ej - E^J (Ej - E^^) 2 (E^ - Ej^) {E^ - E^^y 



-"11 J 

-\E~it 



2 [E^ - Ej^) {E^ - E^^) 



\gr\"\gr\\,.,.5. 



(224) 



AQ{ncccc;t^e,t^e) = 



e~'-®7i* 



+ 



{-itf 



(Ej - E^J [E^^ - Ey) 

-iE.t 



e 1 



[Ery - Ey) (E'-yi - Ey) 



111 

gl 



gV'gr ^77' 



(225) 



AQ{cccnn,kkkc;t^e,t^e) = ^ 



-lE^t 



{E^ — -^71) (-^7 — -^7') 



+ 



^-iE.,,t 



2 (-E-y — -^7') (-^71 — E-yl) 



77 

5i 



„77i „7i7 



(226) 



AQ{cccnn,kkkn;t^e,t^e) = 



71.72 



(-it)'e-'^-'*|g7'^M'g?"^gr'^'^77'^7i72^7i7' ^ 
(Ej - E^^ ) (Ej - Ej^ ) {E^ - E^i ) 



(227) 



45 



Ae{ccncn,kkc;t^e,t^e) = 



c 



-iE-,t 



71,72 L ^ (i;^ - (i;^ - Sy) 



1^771 1 1^7172 1 r^^^^S^^, 



(228) 



AQ{ccncn,kkn;t^e,t^e) = ^ 



A6(ccnnc, fcfccfc;f^e, t^e) = 



(-U)^ e-'^^i* |g7^^|^ \gr'tgV'9r'ri,,.^hYV,. 
2 (£;^ - i;^J [E^^ - E^^) {E^^ - Ey) 



77 

5/ 



771 7i7 

9i9i 



Aq {cnccn, kc;t^e,t^e) = ^ 



77 



77i 7i7 

9i 9i 



2 (-^7 — -^7i) (-^7 — -^7') 



(229) 



(230) 



(231) 



^6 {cnccn, kn; t^e, t^e) = 



e-ii^.* IgJ'^^'' \9rf9r'9r'v,,'V,.,.Vr 



17' 



^6(c?^cnc, fccfc; i e, i e) 



{E^ - E^J {E^ - Ej^) {Ej - E^') 
{-itf 



E 

71^72 



Q—iEyt 



3! (E;^ - Ej^f {E^ - E^^ ) 



{-\tf 



Q—iE^t 



2 {Er^ - E^^f {E^ - E^^f 



{-ity 



Q—iEyt 



2 {E^ - Ej^f {E^ - Ej^ ) 

Iffl'*'^! ^7172 •^77'- 



(232) 



(233) 



Aq {cncnc, knk; t^e, t^e) 



E 

71,72,73 



{-itf 



{-it)' 



3! {E^ - Ey^) {E^ - E^^) {E^ - E^^) 

Q—iE^t 



2 


{E, 


- Ej^) {E^ - 


- ^72) (-^7 ~ 


- E-y^f 


{-it)' 




e~ 






2 


{E^ 


- E^J {E^ - 


- E^^f {Ey 


- Eji) 


{-it)' 




e" 


-lE^t 




2 


{E, 


- EyJ^ (^7 


- E^^) {E^ 


- ^73) 



I ^7172^77173 ^7273 '^77'' 



(234) 



AG{cnncc,kckk;t'e,t^e) = ^ 



(-it) 



2 c 



77 



771 71 7 

9i 9i 



2 (-^7 — -^7')^ (-^71 — -^7')^ 



(235) 



AQ{ncccn,c;t e,t e 



(-it)^e-'^^i'|g7'^M'|gr"1%772'^7 

/ y o / 771 ^2 ^ r:i 771 ^2 ' 



-^77 



71,72 



2 (E7 - E^J^ (E^i - E^J^ 



(236) 



46 




47 



Aq {cncnn, kkkc, kck; t^e, t^e) = 



(-it) 



iE^t U77i|2 



2 e"'-^-!'* 



77 



^771^717 



2 (-^7 — -^71 ) (-^7 — ■^'y' ) 



(248) 



AQ{cncnn,kkkc,knk;t e,t e) 



77 



gj77 ^7 7 rij^^j^Vlil' 



VI .-y, ^ (-^7 - ^71 ) (-^7 - -^72 ) (-E^7 " -^7' )^ 



71)72 



(249) 



AQ{cncnn,kkkn,kck;t'^e,t^e) = 



(-it)-^c IffrT Igj'l'gi "tfi'" V-VVin2Vo2- 
2 (E^ - E^^ ) (E^ - ) (£^ - iJ^/ ) 



(250) 



Ae{cncnn,kkkn,knk,kc;t^e,t^e) = 



77 



772 727 „ „ , 
i/l i/1 '/7i72'/7l7 



(251) 



Ae{cncnn, kkkn, knk, kn; t^e, t^e) 

= E 



(_if)2 e-iiJ^t 1^771 |^|^^77. 

I ffl 5l ^77''77l72''77l73''77l7''77273^7: 



'27 



(-^7 - -E71) (-E7 - -E72) (■E'7 - -E73) (£^7 - -E7') 



(252) 



AQ{cnncn,kckk,kkc;t^e,t^e) = 



(_it)2e-ii^.*|577l|2 1^71 72 |2 1^772 



77 



7i>72 



(-^7 ~ -^71)^ (-^7 ~ ^12^ 



(253) 



Ae{cnncn,knkk,kkc;t\t'e) = J2 '^'^''^ |gf ^1%773 ^7i72^7i73'^77- 

71,72,73 (-^7 - -^71) (^7 - -^72) (-^7 ~ -^73) 



(254) 



AQ{cnnnc,kcck;t^e,t^e) = 

71 



2 (-E'7 — -E71) (-^7 — -^7')^ 



+ 



(-it)2 



2 (-£/7 — Ej' ) (-^71 — -f/'y' ) 



57" 



^771^717 



(255) 



(-it) 



2 c 





2 




2 











A(,{cnnnc,kcnk]t^e,t^e) = — ,r 2 

71,72 (-^7 - -^7') (-^71 ~ -^7') (^72 - ^7') 



(256) 



(_ii)2e-'^^M57"^' 



77 



71,72 ^ (-^"^ ~ -^Ti) (-^7 - -^72) {^12 - Ej'f 



(257) 



Aq {nccnn, kkkc, ck; t^e, t^e) = 

71 



-it) 



2 e^'^v* 



717 

57 



„77 „7 7i„7i7„77 

5i yi 5i yi 



2 (-^7 — -^7') (-^71 — -^7') 



(258) 



Ae{nccnn,kkkc,nk;t^e,t^e) = 



(-it) 



2 e 



717 
57 



71,72 



(-^7 ~ -^7') (-^71 ~ -^7') (-^72 - -^7') 



(259) 



48 



Ae{ncncn,ckc;t^e,t^e) = 



i-it) 



2 e 



'\9T 



111 

9i 



{Ej - Ej^) {Ej^ - Ey) 



(260) 



AQ{ncncn, ckn; t e,t e 



2„ ,3„^ _ (-ip^ e-^^^i* \gr^f Ig^'f Q^gr' V,-y.V-yyV,.-y' 



y1~^, ^ (-^7 ~ ^71 )^ (-2^71 - -^72) (-^71 - -2^7') 



71,72 



(261) 



AQ(ncncn,nkc;t^e,t^e) = — 



(-it) 



2 e~'-^'^i* 



^717 



1^7172 I g'Y^'^ gj^'^ ?777''?772 ^77'%27' 



(-E^7 ~ -^71) (-^71 ~ -^72) (-2^71 ~ -^7') 



(262) 



AQ{ncncn,nkn,c;t e,t e) 



e-i^^-a* 1,97^/^1^ \gJ^->-f l.97^"1%77.^773^7273'577' 



^ 2 

71,72,73 



(-^7 ~ -^71) (-^71 ~ -^72) (-^71 ~ -^73) 



(263) 



AQ{ncncn, nkn, n; t^e, t^e) 

I ffl 9l ^772'7773 ^77''?7273'?727'^73 



= - E 



2 (i;^, - E^J {Ej^ - E^^) [E^^ - Ej^) {Ej^ - E^>) 



(264) 



Ae{ncnnc,kkck,ckk;t^e,t^e) = 



(265) 



AQ{ncnnc,kknk,ckk;t e,t e 



^ (-it)^ e-'^^* Igr^f \grp' I ffZ^'^l' ^772^7173^7273'^- 

1,^2,^3 ^ (-E'7--^7i) (-E'7 ~ ^72) (-E'7 ~ -^73) 



77 



(266) 



^6(?^?^ccn, cc; f^e, f^e) = 



71 



(-it)^ e-'^-'* \gr' f gT' gJ^^gT'^gV' 

2 {E^ - E^, f {E^ - Eyf 



(267) 



A6(nnccn, cn; t^e, f'^e) = 



(_it)2 ^-iE^t |g772 |2 g7^^g7^72g727g77'^^^^,^^^^, ^ 

2 (£;^ - iJ^j ) [Ej - E^^)^ {Ej - Ey ) 



(268) 



AQ{nncnc,ckkk,kck;t'^e,t^e) = 



(-it) 



2 e-'^v* 



77 

gV 



7i7 

5i 



^771^717 



2 (-^7 — E-y')'^ {E-yi — -^7')^ 



(269) 



Aq {nncnc, ckkk, knk; t^e, t^e) = 

71,72 



i-it) 



2C '^V* 



77 

gr 



gr 



2 (i;^ - Eyf {E^^ - Ey) [E^^ - Ey) 



(270) 



Ae{nncnc,nkkk,kck;t e,t e) 



2p ^3 s ^ V- 



(-it) 



2 e-'^v* 



^7.7 



^727 



771 71 7 

5i 5i '7772'?7i 72^777' 



71,72 ^ (-^7 ~ -^7') (-^71 ~ Eyf {E^^ - Ey) 



(271) 



AG{nncnc,nkkk,knk,ck;t^e,t^e) = 



(-it) 



2 e '^v* 



77 



.91 



771 717 

5l 5l '7772^7172 



71,72 



2 (£;^ - ii^7')^ (£"71 - E^') {Ej2 - Ey) 



(272) 



49 



AQ{nncnc, nkkk, knk, nk; t'^e, t^e) 



E 



Ae{nnncc,cckk;t^e,t"^e) — 





2 




2 











{Ej - Ey) {Ej^ - E^) {E^^ - Ey) {Ej^ - Ery>) 

2 (-B-y — -£'71 ) (-^7 — -^7')'^ 



+ 



57" 



„77 „7 71 „7l7„77 

ffi ffi Si Si 



(273) 



(274) 



Ae{nnncc,cnkk]t e,t e) 



E 

71.72 



sr 



„77l„7l72 72 7„77 „ „ 

yi Si Si yi Vjij'Vj2i' 



VI ,7, ^ (^7 ~ -^71 ) (-^7 - ^72 ) (-£^7 ~ ^7' ) 



(275) 



AQ{nnncc,nckk;t^e,t'^e) = 



(-it) 



2 C 



-IE,, I 



Hi' 



Hi Hi Hi 'Hi' ih'iihri 



71,72 



2 (25^ - Ey) [Ey^ - Ej>) {Ej2 — Ey) 



(276) 



AQ{nnncc,nnkk,c;t^e,t^e) = 



(_if)2 e-i^^.* l57''1'57"^9r"^9r"^97'^\772^7l73'^7r 



71 ,72 ,73 



{Ej - £-yJ (Sy, - Ej^) {E^ - E^^) 



(277) 



AQ[cnnnn,kccc;t^e,t^e) = 



71 



(278) 



Ae{cnnnn,kccn;t\th) = -4^' 2 VjiYVy.Y ^ 

2 {Ey ~ E^J {E^ - E^^) {Ey - Ey) 



(279) 



AQ{cnnnn,kncc;t'^e,t^e) = 



71,72 



2 {Ej - Ej^ ) (i;^ - E^^) {E^ - Eyf 



(280) 



AQ{cnnnn,kcnn,kkc;t^e,t^e) = 



(_if)2 e-i^.* IgJ-r^ 1^ ff7"^g7^"^g?"^g7''"r?7.73'^77' 

2 {Ej - Ery^ ) {Ery " E^^ ) {E ry " E.y^ ) 



(281) 



^6 {cnnnn, kncn, kc; t^e, t^e) = 



97" 



„77i„7i72 727 „77 „ ,„ , 

yi yi yi yi n^i-y'^i^i' 



{Ej - Ej^ ) {Ej - E^^ ) {Ej - Ey 



(282) 



^6(crannn, kncn, kn: t^o, f^c) 

( — i^) e ' ^ |ff7" I Hi' 9l ' ffl' "ffl " ^7i72^7l73^7l7'''7727'^737' 



= E 



71 ,72 ,73 



{Ef - -E71) (-E7 - -E72) (-E7 - £^73) (-E7 - £^7') 



(283) 



('-i^'\2 p-i-E^t I 77l|2 772 7273 7371 „7l7„ „ S 

Ae{cnnnn,knnn,kcc;t\t'e) = ^ ' ''^^ , ''^^ ''^^ ''^^ ^773^717 2^7' 

2 (£;^ - ii^^j ) {Ej - E^2 ) {Ej - E^^ ) 



71,72,73 



(284) 



50 



AQ{cnnnn, knnn, knc; t^e, t^e) 



= E 



71 ,72 ,73 ,74 



(-^7 ~ -^71) (-^7 ~ -^72) i^t - -^73) (-£'7 ~ -^74) 



(285) 



Ae{ncnnn,kkcc,ckk;t e,t e) 



2„ +3„\ _ 



(-it) 



2 e^'-^-y'* 



57" 



^717 



^771^717 



2 (£;^ _ _E'^,)2 (_e;^^ _ E^f 



(286) 



AQ{ncnnn,kkcc,nkk;t^e,t^e) = 



(-ii) 



2 e^^^T'* 



^717 



2 / f / 

9l g'J^'^'^ gj'^^ '?77i?7772 



(287) 



^2 e-'^V* 



^6 {ncnnn, kknc, ckk; t e, t e) = 



77 



7i7 

5i 



772 727 „ „ 
yi i/l '/77l"7l72 



^ Hi): ^ 

7, ,79 ^ (-^7 ~ ^1') (-^71 ' -^7') (-^72 - ^7') 



(288) 



AQ{ncnnn,kknc,nkk,ck;t'^e,t^e) = 



(-ii) 



2 6 '^V* 



^7.7 



77 7 72 727 77 

.vi .yi yi yi 7771 '/71 72 



7ii72 



2 (£;^ - £^/) (i;^^ - Ey) [E^^ - Ey) 



(289) 



AQ(ncnnn, kknc, nkk, nk; t'^e, t^e) 

2 



E 

71,72,73 



(-it)2e-'^-'* 



gTI gJ gJ gJ '7771 ^772*7773^7172^7173 



2 {E^ - Ey) [E^^ - Ey) {E^^ - Ey) {E^^ - E^) 



(290) 



AQ[nncnn, ckkc, kck; t^e, t^e 



y (-it)^e-'^^* |g77i|2 1^772 |2 1^7172 |2^^^, 



(291) 



Ae{nncnn,ckkc,knk;t^e,t^e) = ^ 



(-it)2 e-'^^* Igrfgr^gr^gr^gri-r.-fsS^ 



77 



71,72,73 



2 (£^ - Ej^ ) (i;^, - E^^ ) (i;^ - i;T,3 ) 



(292) 



Ae{nnncn,cckk,kkc;t^e,t^e) = 



(_i^^2 e-'^^* K^^r 



77 



771 71 7 

5i 5i 



71 



2 {Ej - E^^f {E^ - Eyf 



(293) 



Ae{nnncn,cckk,kkn;t^e,t^e) = 



(_ii)2 e-ii^.* \grfgr'gr^grgr'v,.,'V,.,' 



71 ,72 2 (-^7 - ^71 )^ (-^7 - -^72 ) {Ej - Ey ) 



(294) 



AQ{nnncn,cnkk,kkc;t^e,t^e) = 



(_if)2e '^-'Ig/ 



iE-,t I 772,2 



77 



771 7l7 
5l 5l '77172*7727' 



71 ,72 2 (-E'7 - -^71 ) (-^7 ~ -^72 ) (^7 - ^I'f 



(295) 



A^{nnncn,cnkk,kkn,kc;t^e,t^e) = 



7i>72 



(_i^)2 e-ii^.* |g7"l'g7" g7'"^gr"g7"'^7i72'/727' 
2 (i*;^ - ) {Ej - E^^ ) (iJ^ - E^> ) 



(296) 



AQ{nnncn, cnkk, kkn, kn; t e, t e 



71 >72 ,73 



( it)^ e Iffl^^l gj^ ffl gj g'P '^7l73^7l7'^7273^727'%37' 

2 {E-y - EjJ {E^ - E^^) {E^ - E^^) {E^ - Ey) 

2 



AQ{nnnnc,ccck;t^e,t^e) = 



i-it) 



2 e T 



7l7 

ffl 



„77 „7 71 „7l7 „77 
i/l i/l i/l i/l 



71 



AQ{nnnnc,ccnk;t^e,t^e) = 



AQ{nnnnc,ncck;t^e,t^e) — 



(-ii) 



2 e"^^i'* 



— E^f ) {E^^ — i?^' ) 



9i 9i 9i 9\ 9i '7772 '/7i72 



{E^ — Ej>) {E^^ - Ey) [E^^ — Ey) 

2 



2 e"'^^'* 



9r 



gV gJ^'gr'-grv^^.v^-y. 



{E^ - Ey) (E-yi - Ey) {Ery^ - Ey) 
2 



Aq {nnnnc, ncnk, ck; t^e, t^e) = 



i-it) 



2 e"'^^'* 



77 



^77 g'J g'^ g'J '7771 ^7772 



7i>72 
2„ +3, 



AQ{nnnnc,ncnk,nk;t e,t e 



E 

71 ,72 ,73 



i-it) 



2 



^737 



(-E'7 ~ -^7') (-^71 ~ -^7') (-^72 - E-y') 



gTi gj ^771 '7772^773 ^7173 ^7273 



Aq {nnnnc, nnck, c; t^e, t^e) = 

■i>7 

A(,{nnnnc,nnnk,cck;t^e,t^e) = ^ 



[E^ — Ey) [E^^ — Ey) {Ej^ - Ey) [E^^ — Ey) 

(_it)2 e-i^.* \gr\' gr'gr^gr'grv,.,.5,y 



71,72,73 ^ (-^7 ~ -^71) (-^7 ~ ^72) (-^7 --^73) 

e-'^-.* |.97"1^97"^97^"^9?"^5?"^772 ^7173^77' 



7ii72,73 



(£■7 - i;7i ) (£'7 - £72 ) (£7 - £73 ) 



AQ{nnnnc, nnnk, cnk; t^e, t^e) 

(_if)2 e-ii^.* 157^*1^ 

.9l ffl ffl ffl '7772'?773'77l 73*771 74 ^72 74*77374 "77' 



= E 



7i,72i73,74 



-l2 ^-iB^t U774|2 771 7172 ^7273 „737^ 
1 I ffl ffl ffl ffl ' 
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